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preparation  and  b)  a  whole-cell  patch  recording  protocol  of  direction- 
ally  selective  ganglion  cells,  the  combination  of  which  is  well- suited  for 
characterizing  complicated  retinal  receptive  fields.  I  will  show  evidence 
that  when  the  inhibitory  synaptic  input  is  markedly  reduced  the  ex¬ 
citatory  input  to  directionally  selective  ganglion  cells  in  turtle  is  itself 
directionally  selective.  This  reduction  was  accomplished  by  a)  voltage 
clamping  to  the  reversal  potential  of  inhibitory  currents  and  b)  by  re¬ 
moving  ATP  from  the  electrodes  which,  in  turn,  blocks  the  inhibitory 
input  over  time. 

This  finding  implies  that  the  necessary  and  sufRcient  conditions  for 
the  computation  of  directional  selectivity  occurs  pre-synaptic  to  the 
ganglion  cell.  For  the  theoretical  side  of  the  Thesis,  and  consistent 
with  the  above  constraint,  I  present  a  model  for  the  direction  selec¬ 
tive  circuit  that  relies  on  directionally  selective  distal  outputs  from 
functionally  independent  oriented  amacrine  cell  dendrites.  I  postulate 
that  these  dendrite  cables  receive  symmetric,  evenly  distributed,  but 
temporally  distinct  excitatory  and  inhibitory  inputs  along  their  length. 
Performance  of  the  model  (tested  with  a  biophysically  and  morphome- 
tricly  detailed  neuron  simulator  program,  SURF-HIPPO)  is  consistent 
with  the  data,  and  generates  testable  predictions. 

I  shall  also  discuss  how  this  model  may  work  in  a  developmental 
context  and,  finally,  implications  for  mote  general  non-linear  spatio- 
temporal  filtering  within  dendritic  trees. 
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Abstract 


The  retina  is  a  good  candidate  for  exploring  the  relationship  between 
neural  computation  and  circuit,  in  particular  given  its  physically  pe¬ 
ripheral  location  and  its  physiologically  central  status.  One  example  of 
a  spatial-temporal  computation  in  the  retina  is  directional  selectivity. 
This  computation  may  rely  on  interactions  within  the  dendritic  tree 
which  are  incrementally  more  complex  than  the  basic  “point  integra¬ 
tion”  and  fire  neuronal  response. 

In  this  Thesis  I  use  experimental  and  theoretical  techniques  to  char¬ 
acterize  the  direction  selectivity  circuit  in  vertebrate  retina.  For  the 
experimental  side  I  have  developed  a)  an  isolated  intact  turtle  retina 
preparation  and  b)  &  whole-cell  patch  recording  protocol  of  direction- 
ally  selective  ganglion  cells,  the  combination  of  which  is  well- suited  for 
characterizing  complicated  retinal  receptive  fields.  I  will  show  evidence 
that  when  the  inhibitory  synaptic  input  is  markedly  reduced  the  ex¬ 
citatory  input  to  directionally  selective  ganglion  cells  in  turtle  is  itself 
directionally  selective.  This  reduction  was  accomplished  by  a)  voltage 
clamping  to  the  reversal  potential  of  inhibitory  currents  and  b)  hy  re¬ 
moving  ATP  from  the  electrodes  which,  in  turn,  blocks  the  inhibitory 
input  over  time. 

This  finding  implies  that  the  necessary  and  sufficient  conditions  for 
the  computation  of  directional  selectivity  occurs  pre-synaptic  to  the 
ganglion  cell.  For  the  theoretical  side  of  the  Thesis,  and  consistent 
with  the  above  constraint,  I  present  a  model  for  the  direction  selec¬ 
tive  circuit  that  relies  on  directionally  selective  distal  outputs  from 
functionally  independent  oriented  amacrine  cell  dendrites.  I  postulate 
that  these  dendrite  cables  receive  symmetric,  evenly  distributed,  but 
temporally  distinct  excitatory  and  inhibitory  inputs  along  their  length. 
Performance  of  the  model  (tested  with  a  biophysically  and  morphome- 
tricly  detailed  neuron  simulator  program,  SURF-HIPPO)  is  consistent 
vdth  the  data,  and  generates  testable  predictions. 

I  sheJl  also  discuss  how  this  model  may  work  in  a  developmental 
context  and,  finally,  implications  for  more  general  non-linear  spatio- 
temporal  filtering  within  dendritic  trees. 
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Chapter  1 

Introduction 


In  this  Thesis  we  investigate  the  generation  of  directionally  selective  re¬ 
sponses  in  the  vertebrate  retina  (Chapter  2).  Our  perspective  will  be  both 
experimental  and  theoretical.  We  shall  use  the  turtle  retina  (and  to  a  lesser 
degree,  rabbit  retina)  as  a  model  system  for  both  vertebrate  retinal  process¬ 
ing  in  particular  and  neural  computation  in  general.  The  major  results  are 
as  follows: 

On  the  experimental  side  (Chapter  3),  we  shall  show  that  at  least  some 
directionally  selective  ganglion  cells  in  the  turtle  retina  receive  directionally 
selective  excitatory  input  ([13]).  This  interpretation  is  supported  by  both 
"oltage  clamt?  recordings  at  holding  potentials  near  inL...  to;  ^  ptic  re¬ 
versal  potentials,  and  evidence  that  inhibitory  input  onto  recorded  cells  is 
attenuated  under  our  recording  conditions.  An  important  implication  of  this 
finding  is  a  constraint  placing  the  location  for  the  computational  elements 
necessary  for  direction  selectivity  onto  cells  pres]maptic  to  ganglion  cells. 

On  the  theoretical  side  (Chapter  4),  we  shall  present  a  model  for  the  pre¬ 
ganglionic  computation  of  direction  selectivity  in  vertebrate  retina  ([12]). 
This  model  is  consistent  with  known  properties  of  the  retina  and  the  afore¬ 
mentioned  results,  and  requires  reasonable  assumptions/constraints  on  what 
is  not  known.  Computer  simulations  of  biophysically  and  morphometricly 
detailed  models  of  amacrine  cells  demonstrate  the  performance  of  the  model, 
and  suggest  that  the  model  is  robiist  for  realistic  variations  in  model  param¬ 
eters. 

This  model  is  based  on  the  sequence-dependent  output  of  a  dendritic 
cable  resulting  from  non-linear  interactions  between  the  cable’s  excitatory 
and  inhibitory  synaptic  inputs.  We  suggest  that  the  biophysical  substrate 
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for  this  model  might  be  found  in  extra-retinal  systems  (Chapter  5),  and 
we  discuss  the  possibility  of  such  computations  (non-linear  spatial-tempor2J 
input  filtering)  in  other  contexts.  We  also  discuss  advantages  that  this 
particular  model  has  for  solving  the  developmental  problem  of  the  breaking 
of  symmetry,  and  what  this  model  suggests  about  ganglion  cell  function  in 
general. 

Finally,  this  project  included  the  development  of  two  new  experimental 
protocols  -  isolation  and  recording  &om  intact  turtle  retina  ([11]),  and  whole 
cell  patch  recording  from  intact  turtle  retina([13]).  These  techniques  are  dis¬ 
cussed  in  Appendix  A  and  6.  This  project  also  included  the  development  of 
a  circuit  simulator  program,  SURF-HIPPO,  that  facilitated  our  theoretical 
investigations.  SURF-HIPPO  is  discussed  in  Appendix  C. 
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Chapter  2 


Directional  Selectivity  in  the 
Visual  System 


2.1  A  Definition  of  Directional  Selectivity 

What  is  directional  selectivity  (DS)  in  the  visual  system?  In  this  thesis  the 
definition  is:  The  response  of  some  visual  pathway  to  a  moving  stimulus  is 
DS  when  some  linear  measurement  taken  over  the  entire  response  (e.g.  the 
integral)  is  significantly  asymmetric  with  respect  to  stimulus  direction.  An 
example,  and  perhaps  t*- '  classical  definition,  would  be  that  a  cell  is  DS  if 
it  consistently  fires  more  spikes  for  a  specific  (preferred)  stimulus  trajectory 
as  con^)ared  to  the  (null)  stimulus  trajectory  that  differs  only  in  sign.*  On 
the  other  hand,  for  a  non-spiking  cell,  an  EPSP  with  a  significantly  larger 
area  for  a  specific  motion  direction  is  defined  as  DS.  A  common  thread  for 
our  definition  is  that  a  DS  system  requires  the  breaking  of  symmetry  over 
the  receptive  field  of  that  system. 

Behaviorly,  of  course,  the  most  important  DS  signals  are  carried  by  spikes 
(we  leave  aside  the  possibility  of  subthreshold  signaling  within  a  subsystem, 
e.g.  an  intra-retinal  DS  channel).  Certainly  the  motivation  for  studying 
DS  starts  with  spike-based  systems,  but  in  this  Thesis  we  are  ultimately 
interested  in  dissecting  the  mechanisms  underlying  a  specific  computation. 
Strictly  speaking,  this  may  be  done  independent  of  considering  the  function¬ 
ality  of  the  specific  computation.  In  the  case  of  retinal  DS,  this  means  that 
all  DS  somatic  ganglion  cell  responses  are  relevant,  since  it  is  straightforward 

'Preferred  or  null  msj  refer  to  either  the  stimulua  or  the  eyitem/ceQ  reeponee,  m 
approprinte. 
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to  apply  known  mechanisms  to,  for  example,  DS  EPSPs  to  yield  the  classi¬ 
cal  signals.  Thus  in  this  Thesis  we  may  consider  motion  asymmetries  which 
may  either  account  directly  for  the  classical  DS  responses,  or,  in  principle, 
be  used  by  the  brain  for  directional  information. 

2.1.1  Directional  Selectivity  versus  Difference 

In  fact,  for  our  analysis  of  DS  we  will  also  consider  a  broader  directional 
difference  (DD)  response  distinction,  one  which  also  requires  a  fimdamental 
asymmetry:  a  system  or  cell  output  is  DD  if  there  is  any  consistent  pre- 
ferred/null  (P/N)  difference  that  is  independent  of  any  relative  time  shift 
of  the  two  waveforms.  For  example,  P/N  output  waveforms  with  equivalent 
averages  but  different  shapes  or  amplitudes  would  constitute  a  DD,  but  not 
DS,  response.  As  described  by  Poggio  and  Reichardt  ([55]),  DS  responses 
also  require  a  nonlinearity  in  the  circuit. 

2.1.2  Consistency  of  Directional  Properties  Over  the  Recep¬ 
tive  Field 

In  either  case,  in  this  Thesis  consistency  means  that  a  single  P/N  axis  is 
dominant  over  the  entire  receptive  field  of  the  cell.  This  constraint  is  mo¬ 
tivated  by  both  the  classical  work  on  DS  (e.g.  [8]),  and  by  the  notion  of 
breaking  of  symmetry  mentioned  earlier.  Recently,  center/surroimd  DS  an¬ 
tagonist’ .  receptive  fields  have  been  reported  in  turtle  retina  ([22]),  but  in 
principle  these  cells  would  be  considered  DS  since  they  still  break  symme¬ 
try.  As  a  counter-example,  a  cell  with  a  classical  excitatory  center/inhibitory 
surround  may  give  a  centrifugal  “DS”  response  for  movement  across  the  re¬ 
ceptive  field  edge  (see  Section  3.1).  However,  since  the  observed  P/N  axis 
rotates  as  a  function  of  stimulus  location,  the  cell  would  not  be  classified  as 
either  DS  or  DD. 

2.1.3  Directionally  Selective  Responses  Emerge  at  Different 
Levels  for  Different  Species 

Directionally  selective  visual  responses  emerge  at  different  levels  of  the  visual 
pathway  in  different  species  -  in  general,  complicated  receptive  fields  tend 
to  appear  at  higher  levels  with  higher  species.  For  exanq>le,  directionally 
selective  x  >onses  are  common  in  turtle  optic  nerve  fibers  (ranging  between 
20%  to  40%  [27],  [20]),  whereas  in  primate,  similar  responses  appear  in 
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relative  abundance  first  in  VI  (however,  a  small  percentage  of  primate  retina 
DS  responses  have  been  reported  [57]  (21]). 

2.1.4  Structure  of,  and  Directional  Response  in  the  Retina 

The  ‘‘sensory  membrane  of  the  eyeball”  is  a  good  candidate  for  study  given 
its  physically  pertpheroZ  location  and  its  physiologically  centra/ status  ([24]). 
The  vertebrate  retina  is  organized  in  several  layers  of  cell  bodies  and  their 
interacting  processes.  Signal  flow  is  both  direct  (perpendicular  to  the  image) 
and  lateral  (parallel  to  the  image)  at  all  levels.  Light  is  transduced  at  the 
photoreceptor  layer,  which  outputs  to  bipolar  and  horizontal  ceUs  within 
the  outer  pleziform  layer.  Bipolar  cell  output  impinges  on  the  mesh  of 
amacrine  and  ganglion  cell  dendrites  within  the  inner  pleziform  layer  (IPL). 
Finally,  ganglion  cell  axons  form  the  optic  nerve.  Each  major  cell  t]rpe  in 
the  retina  has  several  subtypes,  classified  either  anatomically  (e.g.  according 
to  dendritic  tree  shape),  neurochemically  (e.g.  cholinergic,  GABAergic),  or 
physiologically  (e.g.  ON/OFF,  DS,  red/green  opponent). 

In  turtle,  motion  dependent  responses  have  been  reported  at  all  retinal 
layers  [23].  The  focus  of  this  Thesis  is  to  analyze  the  circuitry  responsible 
for  DS  responses  at  the  output  layer  -  the  ganglion  cells  -  of  the  vertebrate 
retina,  iising  turtle  retina  as  an  experimental  model  and  turtle  and  rabbit 
retina  as  the  primary  basis  for  a  theoretical  circuit  model. 


2.2  Theoretical  Requirements  for  Directional  Se¬ 
lectivity  and  Difference 

Motion  detection  is  a  computation  on  tpatially  separated  inputs  over  time. 
Detection  of  motion  direction  (DD)  requires  a  spatial  asymmetry  in  the  cir¬ 
cuit.  Finally,  as  mentioned  earlier,  a  non-linearity  is  necessary  to  transform 
a  DD  signal  to  DS.  Note  that  it  is  also  possible  to  have  a  DS  system  in  which 
there  is  no  DD  signal. 

It  is  useful  at  this  point  to  define  two  broad  classes  of  retinal  DS  mod¬ 
els:  ganglionic  models,  where  the  crucial  nonlinear  interaction  occurs  in  the 
ganglion  cell,  and  pre-ganglionic  models,  where  the  interaction  occurs  prior 
to  the  ganglion  cell  (see  review  in  [39];  in  this  paper  post-synaptic  means 
ganglionic,  and  pre-synaptic  means  pre-ganglionic). 

Thus,  the  specific  questions  we  we  trying  to  answer  here  are: 

•  What  is  the  anatomy  and  connectivity  of  the  DS  pathway? 
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•  What  is  the  crucial  nonlinearity  of  the  DS  pathway? 


•  Where  on  the  DS  pathway  is  the  nonlinearity,  e.g.  is  it  pre-ganglionic 
or  ganglionic? 

The  experimental  portion  of  this  Thesis  deals  with  the  last  point  most 
directly,  while  the  first  two  points  are  more  pertinent  to  the  theoretical 
portion. 


2.3  Experimental  Work  on  Retinal  Directional 
Selectivity 

Over  the  past  thirty  years  there  has  been  a  large  body  of  work  investigating 
retinal  DS  with  a  variety  of  preparations.  DS  retinal  output  was  first  de¬ 
scribed  in  amphibian  (&og,  Maturana,  Lettvin,  McCulloch,  and  Pitts  [50]) 
and  eventually  characterized  in  insect  (fly,  e.g.  Hausen  [34]),  reptile  (turtle, 
e.g.  Lipetz  and  HiU  [44]),  bird  (pigeon,  Maturana  [49]),  and  mammal  (e.g. 
rabbit,  Barlow  and  HiU  [7]).  We  shaU  now  overview  some  of  the  key  findings 
in  rabbit  and  turtle  that  are  pertinent  to  this  Thesis. 

2.3.1  Extracellular  Electrophysiology 

The  classical  eztraceUular  rabbit  experiments  of  Barlow  and  Levick  ([8]) 
described  several  phenomena  related  to  the  DS  response.  Using  both  moving 
slits  and  apparent  motion’  protocols  they  showed: 

•  DS  subunits:  SmaU  non-overlapping  stimulus  regions  within  the  slit- 
mapped  receptive  field  which  were  DS  consistent  with  the  response  to 
fuU  field  stimuli. 

•  Inhibitory  Mechanisms:  Stimulus  at  a  given  point  in  the  receptive  field 
inhibited  later  responses  to  stimuli  in  the  ntiU  direction. 

The  first  finding  suggested  that  the  DS  circuit  elements  for  a  given  gan- 
gUon  ceU  were  replicated  many  times  for  that  ceU.  Historicsdly,  the  most 
explored  interpretation  of  the  second  finding  was  that  the  DS  computation 
relied  on  an  asymmetric  inhibitory  pathway.  However,  these  data  cannot 

*  Apparent  motion  stimuli  are  sequences  of  flashed  spots  in  adjacent  locations,  similar 
to  the  li(hts  in  a  movie  marquee. 
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rule  out  a  model  in  which  excitation  is  asymmetric  and  inhibition  is  sym¬ 
metric  (e.g.  £rom  an  antagonistic  center/surround).  Such  a  possibility  is 
supported  by  the  following  finding. 

Apparent  as}rmmetric  inhibition  was  not  the  only  phenomena  observed 
in  DS;  Barlow  and  Levick  ([8])  also  showed  (and  later  Grzywacz  and  Amthor 
[30]): 

•  Facilitatory  Mechanisms:  A  stimulus  at  a  given  point  in  the  receptive 
field  facilitates  later  responses  to  stimuli  in  the  preferred  direction. 

Note  that  null  direction  inhibition  and  preferred  direction  facilitation 
may  be  accounted  for  by  an  a8}nnmetry  in  either  an  excitatory  (more  pre¬ 
cisely,  facilitatory)  and/or  an  inhibitory  pathway. 

Additional  spatial  and  temporal  parameters  for  the  DS  network  may  be 
inferred  by  the  velocity  tuning  and  size  of  DS  receptive  fields  (e.g.  Wyatt 
and  Daw  [69];  Grzywacz  and  Amthor  [30];  Granda  and  Fulbrook  [28]).  For 
instance,  the  minimum  length  (measured  on  the  retinal  surface)  of  the  lateral 
path  in  the  DS  circuit,  near  the  visual  streak,  is  100  to  300  /xm,  as  derived 
from  apparent  motion  protocols. 

Sequence-dependent  responses  may  be  resolved  for  very  short  separa¬ 
tions,  less  than  10  fan.  (Amthor  and  Grzywacz,  unpublished  data).  This 
short  distance  places  a  constraint  vis-4-vis  the  electrotonic  properties  of  the 

pathway:  passive  cable  properties,  specifically  electrotonic  length  c<-n- 
stants  of  hundreds  of  microns,  are  not  well-suited  for  such  fine  distance 
discrimination.  Thus,  this  result  suggests  that  non-linearities  in  the  path 
are  important. 

The  velocity  of  effective  DS  stimuli  ranges  from  approximately  0.01  to 
10  funinu  (0.1  to  lOO^/s  in  the  visual  field).  These  results  add  constraints 
on  the  synaptic  kinetics  in  the  DS  circuit. 

2.3.2  Intracellular  Electrophysiology 

As  mentioned  earlier,  motion  dependent  responses  have  been  recorded  from 
all  cell  types  in  turtle  retina  ([23]).  These  recordings  suggest  that  the  crucial 
elements  for  DS  are  foimd  throughout  the  retina,  but  two  important  caveats 
must  be  kept  in  mind.  First,  the  existence  of  somatic  DS  signals  prior  to  the 
ganglion  cell  layer  does  not  provide,  by  itself,  evidence  that  these  signals  are 
relevant  to  DS  ganglion  cell  properties.  Second,  and  most  important,  what 
constititutes  a  **true"  DS  response  is  open  to  interpretation,  as  we  discussed 
earlier  in  this  chapter  (see  also  Section  3.1) 
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The  classical  paper  of  Marchiafava  [47]  explored  ganglion  cell  input 
sjrnapses  in  turtle,  for  example  by  examining  the  effects  of  injected  current 
on  DS  responses  (Figure  2.1).  This  paper  analyzed  the  synaptic  mecha¬ 
nisms  underlying  antagonistic  center/surround  responses,  and  the  possible 
relevance  of  these  mechanisms  to  directional  selectivity.  He  concluded  that: 

•  Inhibitory  inputs  to  the  ganglion  cell  with  reversal  potential  close  to 
the  resting  potential  contributed  to  the  antagonism  between  stationary 
stimuli  in  the  center  and  2mnulus  of  the  receptive  field. 

•  Moving  spots  activated  inhibitory  inputs  with  kinetics  and  reversal 
potential  similar  to  those  activated  by  stationary  large  spots. 

•  The  difference  in  kinetics  between  inhibitory  and  excitatory  inputs  is 
sufficient  for  the  reduced  null  direction  response  being  accounted  for 
by  a  correlation  between  inhibition  and  excitation. 

In  concordance  with  this  last  point,  Marchiafava  suggested  that  evidence 
of  motion-evoked  inhibition  meant  that  some  part  of  the  inhibitory  input 
did  not  require  diffuse  stimulation,  as  might  be  inferred  by  the  stationary 
center/surround  properties.  Presumably,  this  interpretation  was  based  on 
the  assumption  that  the  motion-evoked  inhibition  seen  for  the  relatively  long 
stimulus  trajectory  shown  in  the  paper  (991  fan.  Figure  8)  was  present  for 
relatively  short  trajectories,  although  this  not  explicitly  demonstrated, 
now^er,  the  ganglion  cell  motion  responses  shown  in  this  paper  may  also 
be  accoimted  for  by  a  DS  excitatory  input  that  is  concurrent  with  a  non-DS 
inhibitory  input.  A  simulation  of  such  a  mechanism  that  reproduces  the 
salient  features  of  Figure  8  in  [47]  is  shown  in  Figure  2.2. 

The  experiments  of  Watanabe  and  Murakami  in  frog  retina  ([67])  also 
looked  at  the  effect  of  intracellular  current  injection  on  DS  ganglion  cell 
responses.  They  showed: 

•  Motion  stimuli-evoked  IPSP’s  in  ON/OFF,  but  not  ON  nor  OFF  DS 
ganglion  cells,  both  with  no  mjected  current,  and  with  depolarizing 
ctirrent. 

They  concluded  that  a  directionally-dependent  interaction  between  ex¬ 
citation  and  inhibition  at  the  ganglion  cell  could  accoimt  for  the  DS  in 
ON/OFF  cells,  but  that  **The  neural  mechanism  responsible  for  the  differ¬ 
ence  in  EPSP  amplitude  in  ON-  and  OFF-type  cells  probably  lies  presynap- 
tic  to  the  ganglion  cells.” 
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2.3.3  Neurochemistry 

The  neurotransmitters  involved  in  the  DS  circuit  has  been  investigated  by 
pharmacological  protocols  in  rabbit  (Caldwell,  Daw  and  Wyatt  [18],  Ariel 
and  Daw  [6])  and  in  turtle  (Ariel  and  Adolph  [5]).  From  this  work  we  can 
conclude  the  following: 

•  Inhibitory  Mechanisms:  Blockage  of  GABAergic  pathways  reduces  or 
eliminates  DS. 

As  before,  this  result  would  seem  to  support  the  class  of  DS  models  in 
which  inhibition  is  asymmetric.  But  as  we  shall  see  in  Chapter  4,  this  data 
is  also  consistent  with  models  in  which  the  inhibitory  pathway  is  symmet' 
ric.  Also,  recent  results  (Smith,  Grzywacz,  and  Borg-Graham  ([59])  with 
GABAergic  antagonists  in  turtle  show: 

•  Inhibitory  Mechanisms:  For  about  50%  of  all  DS  cells,  DS  is  main¬ 
tained  or  reversed  when  GABAerpc  pathways  are  blocked. 

These  experiments  are  demonstrated  in  Figures  2.3  and  2.4.  This  result 
is  similar  to  that  reported  previously  in  fly  (BulthoiF  and  BiilthoiT,  [16]).  As 
we  discuss  later,  existence  of  DS  without  inhibition  places  strong  constraints 
on  the  DS  circuit  model. 

2.3.4  Aratc 

Physiologically  identified  DS  ganglion  cells  have  been  stained  in  rabbit  (Amthor, 
Oyster,  and  Takahashi  [4])  and  turtle  (Jensen  and  DeVoe  [38]).  A  clear  re¬ 
sult  of  this  work  is  that: 

•  DS  Morphology:  The  dendritic  trees  of  DS  ganglion  cells  are  not 
aligned  with  their  P/N  axes. 

Thus  the  morphometric  substrate  for  DS  is  not  immediately  obvious 
from  the  histology. 

2.4  Theoretical  Work  on  Retinal  Directional  Se¬ 
lectivity 

Inspired  in  part  by  the  correlation  models  of  Hassenstein  and  Reichardt 
([33])  for  motion  detection  in  fly,  several  models  for  retinal  DS  have  been 
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proposed  which  describe  the  location  and  biophysical  mechanisms  for  the 
necessary  (Section  2.2)  spatial  asymmetry  and  time-dependent  nonlinear 
interaction  (Figure  2.5). 

Barlow  and  Levick  ([8])  considered  both  asymmetric  lateral  inhibitory 
and  excitatory  pathways  in  the  outer  plexiform  layer,  with  a  nonlinear  inter¬ 
action  at  bipolar  ceUs  between  these  pathways  and  direct  excitatory  path¬ 
ways.  Others,  including  Torre  and  Poggio  ([63])  and  Koch,  Poggio  and  Torre 
([39]),  suggested  that  the  lateral  pathway  might  be  mediated  by  amacrine 
cells,  among  other  possibilities.  They  showed  that  the  interaction  between 
an  asymmetric  lateral  synaptic  inhibition  and  symmetric  synaptic  excita¬ 
tion,  possibly  on  the  ganglion  cell  membrane  itself,  could  provide  the  nec¬ 
essary  nonlinearity  for  DS.  In  particular,  Torre  and  Poggio  showed  how  the 
shtmting  action  of  inhibitory  synaptic  conductances  interacts  non-linearly 
with  other  inputs.  These  models  do  not  explicitly  define  the  mechanism  of 
the  delay,  other  than  to  point  out  that  a  mechanism  which  has  a  low-pass 
filter  characteristic  (as  had  been  suggested  in  fly  DS  literature),  or  slower 
inhibitory  synaptic  kinetics,  might  suffice. 

We  note  that  in  these  circuit  architectures  the  only  directional  signal 
available  is  strictly  DS,  under  any  circumstances.  This  is  because  the  in¬ 
teraction  between  the  asymmetric  and  symmetric  pathways  is  immediately 
nonlinear.  If  the  nonlinearity  is  blocked  (e.g.  blocking  inhibition  for  the 
inhibitory  model),  then  no  interaction  can  take  pl£M;e:  there  will  be  no  DD 
signal  at  a3.'  Thi-*  directly  relevant  to  the  fin!:’.L..ji  of  DS  with 

blocked  inhibition  that  was  mentioned  in  Section  2.3.3. 

Koch,  Poggio  and  Torre  ([41])  (also  O’Donnell,  Koch  and  Poggio  [53]) 
examined  the  electrotonic  structure  of  a  cat  ganglion  cell  and  a  putative  DS 
rabbit  ganglion  cell  ([4]),  respectively,  in  detail.  They  showed  that  the  den¬ 
dritic  tree  of  the  ganglion  cell  was  well  suited  for  local  interactions  within  the 
tree  between  an  excitatory  input  and  an  inhibitory  input  that  has  a  strong 
shunting  component  (large  conductance  change  with  synaptic  reversal  po¬ 
tential  near  rest).  The  conclusion  was  that  the  computational  substrate  for 
subunit  response  was  possible  within  the  tree,  e.g.  supporting  a  ganglionic 
model. 

2.4.1  Implications  of  Ganglionic  Directionally  Selective  In¬ 
teractions 

Although  not  explored  in  the  work  cited  above,  a  ganglionic  DS  model  also 
predicts  the  arrangement  of  DS  subunits  with  respect  to  the  ganglion  cell 
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dendritic  tree,  and  the  relative  response  to  preferred  direction  motion  which 
traverses  the  entire  receptive  held,  given  certain  constraints.  Assuming  that 
DS  interactions  occurr  at  random  locations  in  a  radially  symmetric  ganglion 
cell  dendritic  tree,  then  preferred  direction  inhibition  is  predicted,  despite 
the  locality  of  the  interactions  (see  Figure  2.6).  For  long  distance  motion, 
this  arises  because  the  distal  tree  on  the  preferred  direction  side  is  activated 
last.  In  order  to  get  to  the  soma,  the  excitation  from  this  area  must  traverse 
dendritic  regions  in  which  the  locally  non-correlated  delayed  inhibition  re¬ 
mains  activated.  This  phenomena  would  be  expected  to  limit  the  spatial 
and  temporal  bandwidth  of  effective  preferred  direction  stimuli.  Likewise, 
subtinits  on  the  preferred  side  of  the  receptive  field  would  be  expected  to 
have  a  smaller  preferred  response  than  that  firom  subunits  on  the  null  side. 

Barlow  and  Levick  ([8])  compared  subunit  responses  to  the  (stationary) 
spot-mapped  receptive  field.  If  we  assume  that  that  the  latter  noapping  is 
congruent  with  the  dendritic  tree,  then  their  results  (Figure  4  in  [8])  do  not 
support  the  prediction  of  preferred  direction  inhibition.  In  this  figure,  for 
example,  the  preferred  direction  response  for  subunits  d  *-*  e,  i  j,  and 
n  ^  o  (all  on  the  preferred  side)  are  similar  to  that  for  o  6,  /  <-►  p, 
and  k  *-*  I  (all  on  the  null  side),  respectively.  On  the  other  hand,  there 
was  a  significsuat  null  direction  response  within  the  three  null  side  subunits 
(anomalous  response,  see  Section  4.3.4).  Further  work  will  be  necessary  to 
determine  whether  these  observations  are  significant. 

Correlation  models  also  precact  v';lor5+-  .oning  as  ..  ■  ction  of  the  dcl?y 

and  length  of  the  asymmetric  pathway,  At  and  /,  and  for  ganglionic  model 
versions  there  is  an  inherent  tradeoff  between  bandwidth  of  null  direction 
motion  and  the  above  preferred  direction  inhibition.  Consider  the  simple 
version  shown  in  Figure  2.5  (B),  where  /  refers  to  the  spacing  between  the 
inhibitory  and  excitatory  pathways.  If  we  assume  that  the  synaptic  temporal 
transfer  functions  are  given  by  S{t)j  then  for  a  moving  spot  the  only  null 
stimulus  is  motion  towards  the  left  with  velocity: 

Vnull  = 

All  other  velocities,  even  towards  the  left,  will  pass.  To  increase  the  rejection 
bandwidth  we  can  add  some  width  r  to  the  inhibition  impulse  response.  Now 
the  rejected  null  direction  velocities  ue  in  the  following  range: 

//(Af  +  t)  <  <  I/At 

However,  since  the  structure  in  which  the  interactions  takes  place  is  symmet- 
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ric,  letting  t  >  0  allows  the  preferred  direction  inhibition  described  above 
(see  Section  4.3.5). 

2.4.2  Starburst  Amacrine  Cell  Models  for  DS 

Recently,  Vaney  and  co-workers  ([66]  [64])  have  expanded  on  ideas  from 
Masland,  Mills,  and  C2issidy  ([48])  and  suggested  a  specific  cell  type  in 
rabbit,  the  starburst  cholinergic/ GAB Aergic  amacrine  cell,  as  providing  an 
asymmetric  excitatory  pathway  for  the  DS  circuit.  They  further  suggested 
that  starburst  dendrites  also  mediate  a  lateral  inhibitory  pathway  in  the  DS 
circuit,  since  these  cells  contain  GABA  and  the  tips  of  adjacent  starburst 
cells  are  in  close  proximity.  In  their  “co-transmission’’  model  (the  same 
cell  releases  both  GABA  and  acetylcholine),  inhibitory  coimectivity  from 
starburst  cells  is  synunetrical,  in  contrast  to  the  asymmetric  excitation.  On 
the  other  hand,  the  locus  and  biophysical  mechanisms  for  the  computation 
of  DS  are  not  specified  in  this  model,  although  they  mention  both  ganglionic 
and  pre-ganglionic  alternatives. 


2.5  Summary 

We  have  described  a  working  definition  of  the  DS  (and  DD)  visual  response 
and  summarized  the  specific  goals  of  this  Thesis.  We  have  also  given  an 
overview  of  both  experimentaJ  i  id  theorel.-iJ  wri-k  t  visual  system  DS, 
especially  as  regards  to  the  vertebrate  retina. 

It  is  important  to  be  precise  with  respect  to  the  visual  system  phenomena 
tmder  scrutiny,  since  quite  different  hypotheses  may  derive  from  subtle  dif¬ 
ferences  in  experimental  or  theoretical  conditions.  For  example,  the  retinal 
DS  model  described  in  Chapter  4  includes  predictions  for  (at  least)  DD  sig¬ 
nals  within  the  circuit  that  underly  the  DS  output.  These  DD  signals  will 
identify  where  in  the  retina  the  specific  computation  of  DS  first  appears. 
Such  a  DD  finding  may  be  used,  in  principle,  to  rule  out  experimentally 
models  which  predict  that  a  strictly  DD  signal  will  not  be  found  elsewhere 
in  the  circuit,  under  any  conditions.  While  the  sloppyness  of  real  neural  cir¬ 
cuitry  would  generate  some  degree  of  asymmetric  responses  in  general,  the 
appropriate  experimental  protocol  should  be  able  to  determine  the  presence 
or  absence  of  a  significant  DD  response  that  is  relevant  to  these  alternate 
models. 
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Fig  I.  The  «;■«<  of  miric  "'•J-f  cuirem  in;«iion  upon  ihe  phoioretponici 'to  g  ma>in|  ipoi.  Ttv: 
trajectory  ol  i...  tpoi  (9SI  jiin  at  }.SM,'mtecl  'axi  along  the  prercrre^null  aaii.  The  retpon^n  to  move- 
mem  in  Ihe  preferred  and  null  direction!  (marked  by  a  icar)  arc  ihows  at  the  left  and  the  nghi  side  of 
hgurc,  respectively.  At  the  top.  the  duration  and  the  direction  of  movement  arc  indicated.  Control 
indicates  responses  recorded  at  testing  potential.  Upper  and  lower  recoids  show  the  photoresponws 
dunng  artif.cial  polariulion  of  the  cell  membrane  obtained  with  I^c  curicnl  intensity  indicated  by  the 
numbers  above  the  records  on  the  left.  Zero  time  coincides  wnh  the  onset  of  the  spot  tnovemenu 


Figtire  2.1:  Protocol  by  Marchiafava  (Figure  8  from  [47])  investigating  the  effect  of 
current  injection  on  DS  responses  in  turtle  retinal  ganglion  cells.  An  important  feature 
of  this  protocol  is  that  for  both  directions  of  the  stimulus,  the  hyperpolarizing  current  is 
sufficient  for  the  synaptic  potentials  to  stay  below  the  resting  potential.  The  small  EPSP 
for  the  null  stimulus  under  hyperpolarization  implies  that  a  null-direction  correlation  of 
inhibition  and  excitation  does  not  underiy  DS  in  this  cell. 
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Figure  2.2:  Simulation  of  lumped  cell  with  OS  excitatory  input  and  non-DS  inhibitory 
input.  Cell  parameters,  includinf  tq,  synaptic  conductances’  kin''*'cs,  sizes,  and 
reversal  p.-  .ntials,  and  vottage-dependent  channels  underlying  the  spike  v.cic  Cnosen 
to  recapitulate  data  from  [47]  reproduced  in  Figure  2.1.  Note  that  when  the  cell  is 
hyperpolarized  below  the  reversal  potential  of  the  inhibitory  synapse  (set  here  at  Brest), 
the  activation  of  both  the  excitatory  and  inhibitory  synapses  in  the  preferred  direction 
combine  to  produce  an  "EPSP"  that  is  higher  than  that  produced  by  the  inhibitory 
synapse  alone  in  the  null  direction,  as  seen  in  Figure  2.1.  Also,  this  choice  for  the 
inhibitory  reversal  potential  means  that  the  membrane  must  be  away  from  Brest  in  order 
to  see  the  activation  of  the  "silent"  inhibition.  In  order  to  account  for  the  distinct  spike 
thresholds  (e.g.  high  threshold  spikes  for  preferred  stimulus  when  cell  is  depolarized), 
two  putative  Na'^  channels  are  included  with  midpoints  of  activation  set  to  -60  mV 
and  -25  mV,  respectively.  These  channels  may  account  for  the  multiple  spike  heights 
reported  in  DS  ganglion  cells  in  turtle  by  Jensen  and  DeVoe  ([38]),  e.g.  differential 
activation  of  more  than  one  type  of  Na^  channel.  Such  diversity  of  Na'^  channels  has 
been  reported  in  other  neurons  and  has  been  predicted  by  modelling  studies  ([9]). 
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Raspon«« 


Figiire  2.3:  Smith  et.  «l.  ([59])  tested  the  sensitivity  of  DS  to  picrotoxin  (a  GABAergic 
antagonist)  in  isolated  intact  turtle  retina,  using  the  extracellular  preparation  described 
in  Appendix  A.  In  some  cells  DS  was  maintained  despite  saturating  levels  of  picro¬ 
toxin  in  the  superfusate  (50  to  200  /iM),  as  demonstrated  in  the  example  cell  above. 
Motion  stimulus:  200/im  X  600/im  bar  at  various  angles,  speed  ranging  from  1  to 
5  |im/ms,  centered  on  the  spot-mapped  receptive  field.  In  this  figure,  contrast  is  de¬ 
fined  as  (/  -  b)/&,  where  /  is  the  foreground  light  intensity,  and  6  is  the  background 
light  intensity. 


26 


Figure  2.4:  Smith  et.  al.  also  reported  that  in  some  cells,  blockage  of  GABAergic 
pathways  resulted  in  a  reversal  of  the  OS  P/N  axis  along  with  the  expected  increase  in 
overall  response.  Experimental  and  stimulus  conditions  similar  to  that  in  Figure  2.3. 


Figure  2.5:  Correlation<type  models  for  the  computation  of  DS  are  typified  by  that 
proposed  for  the  fly  and  vertebrate  retina.  Versions  of  this  model  have  considered 
both  an  all-excitatory  interaction  (A)  or  an  excitatory /inhibitory  interaction  (6)  (e.g. 
asymmetric-veto  or  AND-NOT  models).  In  the  diagram  the  nonlinear  interaction  be¬ 
tween  the  ''irect  and  asymmetric  delayed  inputs  i;.  expressed  Ui.  .  rnultiplic.'''i''n,  but 
this  particular  nonlinearity  is  not  crucial  to  the  model,  in  all  figures  open  and  closed 
circles  are  excitatory  and  inhibitory  inputs,  respectively.  Note  that  a  position-dependent 
delay  is  applied  before  the  nonlinearity.  For  both  (A)  and  (B)  the  preferred  response 
is  for  motion  to  the  right.  For  rightward  motion  in  (A),  the  delayed  lefthand  input  is 
correlated  with  the  un-delayed  righthand  input;  for  leftward  motion  the  delay  amplifies 
the  temporal  separation  at  the  nonlinearity.  For  rightward  motion  in  (B),  the  delayed 
righthand  input  is  correlated  with  the  un-delayed  lefthand  input,  and  the  inversion  of 
the  righthand  input  cancels  out  the  lefthand  input;  for  leftward  motion  again  the  delay 
amplifies  the  temporal  separation  at  the  nonlinearity,  in  this  case  allowing  passage  of 
the  lefthand  input  (assuming  that  the  nonlinearity  is  not  strictly  multiplicative). 
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Time 


Preferred  Direction  Stimulue 


Figure  2.6:  DS  models  in  which  the  ganglion  cell  dendritic  tree  is  the  site  for  DS 
excitatory-inhibitory  interactions  (correlation  of  local  inputs  in  the  null  direction)  predict 
inhibition  of  the  late  phase  of  long  distance  preferred  direction  motion  and  subunit 
preferred  response  that  depends  on  location,  assuming  that  inputs  occur  throughout 
the  tree.  This  inhibition  is  from  previouriy  activated  inhibition  (e.g.  at  C)  that  is  non- 
correlated  with  local  excitation,  but  that  is  nonetheless  both  on-the-path  and  correlated 
for  the  later  excitation  (e.g.  at  D). 
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Chapter  3 

Directionally  Selective 
Ganglion  Cells  in  Turtle: 
Experimental  Results 


We  shtJl  now  present  data  from  whole*cell  patch  (intracellular)  recordings 
of  ganglion  cells  in  turtle  retina.  As  we  saw  in  the  last  chapter,  various 
models  for  retinal  directional  selectivity  have  different  predictions  as  to  the 
properties  of  synaptic  input  onto  directionally  selective  ganglion  cells.  The 
recordings  presented  here  shor^  that  smr  .  mglion  cells  receive  DS  ezc'*-*). 
tory  input:  therefore  the  necessary  elements  for  the  DS  computation  must 
occur  in  cells  pre-synaptic  to  the  ganglion  cell. 

Further,  we  present  on-ceU  patch^  (extracellular)  and  subsequent  whole¬ 
cell  patch  recordings  from  the  same  cells  in  which  extracellularly-determined 
DS  properties  are  coiuistent  with  DS  properties  of  the  excitatory  input. 
These  recordings  suggest  that  the  whole-cell  patch  protocol  does  not  intro¬ 
duce  artifactusd  DS  properties.  DetaUs  on  the  tissue  preparation,  stimulus 
generation,  and  experiment  control  apparatus  are  described  in  Appendix  A. 
Details  of  the  whole-cell  patch  technique,  and  basic  properties  of  cells  recorded 
with  this  technique,  are  described  in  Appendix  B. 

'  On-cell  petch  recordings  ore  nude  with  •  potch  electrode  in  o  tight  seal  configuration 
onto  the  cell,  with  the  cell  membrane  intact. 
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3.1  Stimuli  for  Characterizing  Directional  Selec¬ 
tivity 

A  variety  of  light  stimuli  with  different  motion  components  may  be  used 
to  explore  directional  properties.  Depending  on  the  spatial  and  tempored 
properties  of  the  stimtilus,  each  t]rpe  emphasizes  different  aspects  of  motion 
sensitivities. 

3.1.1  Gratings  Versus  Bars,  Temporal  and  Spatial  Window¬ 
ing 

As  mentioned  in  the  last  chapter,  previous  work  investigating  directional  vi¬ 
sual  responses  have  used  both  aperiodic  stimuli  (e.g.  moving  bars  or  spots, 
apparent  motion  protocols)  and  quasi-periodic  stimuli  (e.g.  spatially  and 
temporally  windowed  gratings).  Both  classes  of  stimuli  have  strong  (direc¬ 
tional)  motion  components.  Aperiodic  stimuli,  however,  emphasize  the  tran¬ 
sient  (in  time)  and  wide  (in  space)  components  of  the  DS  circuit’s  transfer 
fimction,  while  the  quasi-periodic  stimuli  tend  to  emphasize  the  steady-state 
and  narrow  (in  space)  components.^ 

In  principle  any  realizable  DS  circuit  will  show  a  DS  response  to  periodic 
stimuli.  However,  there  are  several  candidates  for  inactivating  or  adaptive 
mechanisms  along  any  neural  pathway  which  tend  to  emphasize  the  tran¬ 
sient  response  of  th'  system,  as  shown  in  Figure  3.1.  Examples  include  ei- 
thor  inactivation  or  adaptation  of  phototransduction,  synaptic  release  (e.g. 
transmitter  depletion),  synaptic  receptors,  and  of  voltage  or  Co'*"*' dependent 
channels  (e.g.  adaptation  of  repetitive  firing,  channel  inactivation  [36]).  In¬ 
hibitory  feedback  may  also  be  important  in  retinal  function  [68].  Also,  in 
retina,  inhibitory  feedforward  pathways  are  relatively  slower  than  their  ex¬ 
citatory  counterparts. 

The  inhibitory  pathways  in  the  retina  are  not  only  relatively  slower  than 
the  parsJlel  excitatory  pathways,  but  they  tend  to  have  wider  spatial  range 
as  weU.  To  the  extent  that  a  DS  circuit  depends  on  inhibitory  input,  these 
properties  would  tend  to  generate  a  more  asymmetric  response  to  stimuli 
with  larger  size  (e.g.  in  terms  os  moving  spot  path  length  or  moving  grating 
aperture). 

’In  this  Thesis,  the  “stesdy-stste”  response  is  defined  similar  to  that  nsed  in  electrical 
engineering  practice,  i.e.  the  response  of  a  system  at  long  times  relative  to  the  start  of  a 
periodic  input. 
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Stimuli  with  a  relatively  smaller  size  have  another  disadvantage  when 
it  comes  to  characterizing  receptive  field  asymmetries.  We  expect  a  DS 
cell  to  have  a  consistent  asymmetry  over  its  receptive  field,  as  we  discussed 
earlier.  The  classical  antagonistic  center-surround  receptive  field,  while  sym¬ 
metric  overaU,  includes  local  asymmetries  (Figure  3.2)  which  could  give  DS 
responses  to  spatially-limited  stimtili  (assuming  non-linear  center- surround 
interactions).  A  harbinger  of  such  a  situation  would  be  a  strong  DS  sensi¬ 
tivity  to  aperture  size  (e.g.  DS  for  apertures  below  some  maximum  size).  A 
consistent  DS  receptive  field  would  not  show  such  a  sensitivity  to  the  spa- 
tial  change  in  a  larger  stimuli  (although  a  concomitant  increase  in  temporal 
characteristics  might  reduce  the  DS  of  the  complete  response).  Such  mech¬ 
anisms  may  explain  some  previous  findings,  e.g.  spatially-limited  DS  cells 
prior  to  ganglion  cells  [23].  An  alternate  explanation  for  the  spatial  lim¬ 
itation  may,  however,  be  center /surround  DS  antagonistic  receptive  fields 
[22]. 

We  have  used  both  gratings  and  bar  stimuli  in  the  experiments  presented 
here.  Although  DS  responses  were  found  using  both  types  of  stimuli,  prelim¬ 
inary  results  suggest  that  the  most  effective  DS  stimuli  under  the  recording 
conditions  described  in  Appendices  A  and  B  are  spots  (e.g.  100  to  2Q0iaa  in 
diameter)  with  a  long  path  (>  SOO^on,  drifting  &om  0.5  to  5  /im/ms).  Long 
moving  bars  (e.g.  1000  x  200/im)  do  not  seem  to  be  as  effective,  possibly 
because  activation  of  the  inhibitory  surround. 

Moving  gracing  stimuli  were  presented  within  a  square  aperture  (100  to400;im 
sides,  oriented  along  the  motion  axis).  Gratings  were  exposed  for  2  second 
“burst”  windows  from  a  dark  background;  depending  on  the  temporal  fre¬ 
quency  of  the  grating,  each  burst  included  1  to  10  cycles  of  the  grating. 

The  spatial  phase  of  the  gratings  was  set  such  that  at  the  start  of  each 
temporal  window  an  ON  edge  appeared  at  one  side  of  the  aperture.  The 
spatial  period  of  aU  gratings  was  at  least  800  pm,  so  the  aperture  was  com¬ 
pletely  dark  at  the  start  of  the  temporal  window.  Responses  from  a  series 
of  windows  at  a  given  orientation  were  then  averaged  together,  maintaining 
zero  relative  phase.  This  protocol  provided  both  a  transient  response  to 
the  first  ON  edge  drifting  across  the  aperture,  and  the  first  cycles  of  the 
“steady-state”  response.  Although  the  rotating  square  aperture  meant  that 
slightly  different  areas  of  the  retina  were  illuminated  depending  on  orien¬ 
tation,  the  illuminated  areas  for  motions  of  opposite  directions  (relevant  to 
the  DS  determination)  was  identical. 

As  described  in  Appendix  A,  100%  contrast  (point  on  monitor  either  full 
on  or  full  off)  stimuli  were  used  in  all  experiments. 
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Figure  3.1:  Mechanisms  within  neural  pathways  include  many  which  are  time- 
dependent.  These  dependencies  may  have  both  activating  and  inactivating 
components,  but  any  inactivating  link  in  a  serial  path  will  shut  ofF  the  signal  in 
the  path,  independent  of  any  activating  links.  The  net  result  will  be  to  attenuate 
the  system  response  to  continued  presentation  of  the  input. 
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Figure  3.2:  The  relatively  slower  kinetics  of  feedforward  inhibitory  pathways  result  in 
local  DS  regions  in  an  overall  symmetric  center /surround  ceil,  assuming  that  the  center- 
surround  interaction  is  non-linear.  The  preferred  direction  of  a  region  is  centripetal  with 
respect  to  the  receptive  field  center,  and  thus  is  dependent  on  the  location  of  the  region. 
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3.2 


Protocols  for  the  Identification  and  Charac¬ 
terization  of  Directionally  Selective  Cells 

3.2.1  Basic  Protocols 

In  each  the  experiments  described  here,  basic  DS  properties  were  determined 
as  follows: 

1.  Single  cell  response  were  obtidned  with  extracellular  (spiking)  on-cell 
patch  recordings  using  flashing  spots  (typically  400  /xm  square)  at  1 
Hz,  centered  either  at  the  electrode  tip  location  or  roughly  midway 
between  the  electrode  tip  and  the  visual  streak,  when  visible). 

2.  Receptive  fleld  spot  map  was  then  obtained  extracellularly  using  (typ¬ 
ically)  100  fan.  square  spots  flashing  at  1  Hz,  stepping  over  a  5  by  5 
square  gr  '  (  00  -  200  ftm  distance  between  grid  points)(protocol  type 
“RF”).  The  center  of  the  resulting  spot  map  was  estimated  by  eye. 

3.  Using  the  center  of  the  receptive  field  as  the  center  of  rotation,  the 

extracellular  responses  to  a  series  of  moving  gratings  or  spots,  moving 
at  either  8  or  16  directions,  were  collected.  A  DS  cell  was  defined  as 
one  in  which  a  chosen  single- valued  response  parameter,  summed  over 
at  most  three  (for  8  direction  experiments)  or  five  (for  16  direction  ex¬ 
periments)  adjacen-:  ':ns  ovc  .  all  trials,  was-significantly  greater 

than  the  summed  response  to  the  complementary  orientations  (Stu¬ 
dent  T  test,  p  <  0.05).  For  example,  for  a  16  direction  experiment, 
the  (preferred)  value 


3 

P=  T,  ^(^0  +  (*  22.5")) 

i=-3 

was  compared  with  the  (null)  value 

3 

N=  ^(^0  +  (*■  X  22.5")  180") 

«=-3 

for  the  16  values  of  Oo,  where  R{9)  is  the  response  for  direction  0.  When 
P  was  statistically  significantly  greater  than  N,  a  direction  index 

DI,^m  =  {P-N)/(P  +  N) 
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was  then  computed  (see  Section  4.3.2).  The  preferred  direction  of 
the  cell,  $p,  was  defined  as  the  value  of  which  gave  the  maximum 
statistically  significant  value  of  The  degrees  of  freedom  for  the 

Student  T  test  was  defined  as  the  number  of  trials  (2  to  8)  multiplied 
by  the  niimber  of  grouped  orientations  (1  to  5).  These  protocols, 
which  also  generate  an  orientation  tuning  for  the  cell,  are  defined  as 
type  “WF”. 

4.  Whole  cell  (intracellular)  access  was  then  obtained  and  WF  protocols 
were  used  to  determine  if  extracellular  DS  properties  (DS  or  not  DS, 
preferred  direction  and  DS  tuning)  were  maintained  over  time,  under 
both  current  and  voltage  clamp  conditions. 

Single- valued  response  parameters  used  for  DS  determination  (A(9))  in¬ 
cluded  spike  count,  integral  of  membrane  potential  minus  Brest  (for  current 
clamp  protocols),  and  integral  of  the  holding  current  minus  the  total  clamp 
current  (for  voltage  clamp  protocols;  note  that  negative  clamp  current  im¬ 
plied  an  excitatory  event). 

Ideally,  extiacelluax  recordings  were  stable  enough  to  perform  the  DS  de¬ 
termination,  but  in  many  cases  either  whole-cell  access  wsis  obtained  prior 
to  completion  of  the  DS  protocol,  or  the  extracellular  recordings  were  too 
noisy  for  evaluation  in  the  limited  time  available.  Likewise,  in  the  ideal 
cose  wl  ole-cell  access  provided  clear  and  l^hle .  spikes  so  that  r  'apitu- 
iated  protocob  (using  spike  count  as  the  response  parameter)  co\ild  be  more 
directly  compared  with  the  previous  extracellular  protocob.  However,  in 
some  celb  spikes  disappeared  within  minutes,  and  either  synaptic  potentiab 
(via  current  clamp)  and/or  synaptic  conductances  (via  voltage  clamp)  were 
evaluated  for  DS  properties.  In  thb  case,  either  the  DS  properties  of  the 
waveforms  were  obvious  to  the  eye,  or  a  similar  student-T  test  as  above  was 
applied  using  the  appropriate  waveform  integral  as  the  dependent  variable. 

3.2.2  Advanced  Protocols 

For  many  identified  DS  cells,  particularly  when  the  recordings  were  suffi¬ 
ciently  stable,  more  specific  DS  measiuements  were  made.  These  included 
the  following: 

•  In  order  to  replicate  the  data  of  Marchiafava  directly  (Figure  2.1)  and 
thus  highlight  ganglion  cell  inhibition,  preferred/nuU  runs  using  spots 
or  gratings  under  various  current  clamp  conditions  were  run  as  soon 
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as  possible  and  after  at  least  30  nainutes  of  whole-cell  access  (pro¬ 
tocol  type  “PN”).  These  included  a  control  run  (current  clamp  at  0 
pA),  and  runs  with  a  constant  depolarizing  or  hyperpolarizing  currents 
(typically  to  about  ±  10  to  30  mV  around  rest,  i,e.  ±  10  to  30  pA). 
Ideally,  if  extracellular  DS  was  determined  then  the  first  set  of  these 
runs  was  done  immediately  alter  whole-cell  access  in  order  to  avoid 
IPSP  attenuation.  Otherwise,  a  WF  protocol  was  done  immediately 
after  whole-cell  access  and  then  (if  the  cell  was  DS)  the  first  run  of 
this  protocol  was  done. 

•  In  order  to  expl  >re  velocity  tuning  and/or  spatial  dependence,  various 
WF  protocols  were  used  with  different  temporal  and  spatial  parame¬ 
ters,  e.g.  various  spot  speeds  or  grating  temporal  frequencies,  or  spot 
path  lengths  or  grating  aperture  sizes. 

•  In  order  to  evaluate  DS  edge  effects  of  symmetrical  receptive  field 
properties,  and  to  measure  DS  subunit  properties,  a  series  of  5  P/N 
runs  with  gratings  or  spots  of  limited  extent  were  used,  with  stimtilus 
centers  spaced  on  typically  100  /xm  centers  along  the  P/N  axis  and  the 
entire  series  centered  on  the  receptive  field  center. 

•  In  order  to  map  DS  subimit  properties  more  fully,  WF  protocols  using 
spatially  limited  moving  spots  or  gratings  were  \ised  over  a  3  by  3 
grid  (100  -  200  fim  distance  between  grid  poijavs),''with  .atire  grid 
centered  on  the  spot- mapped  receptive  field  center. 

•  In  some  cases,  RF  protocols  were  repeated  at  various  intervals  to  check 
the  stability  of  the  spot-mapped  receptive  field. 

For  some  cells,  simple  size  dependence  for  flashing  spots  was  determined 
using  50  to  1600  ftm.  square  spots,  centered  on  the  spot-mapped  receptive 
field.  This  protocol,  called  type  ‘*SP’’,  was  especially  useful  for  observing 
the  evolution  of  inhibitory  post  synaptic  potentials  (IPSPs)  during  whole¬ 
cell  access.  Because  of  its  importance  in  the  interpretation  of  intracellular 
DS  recordings,  we  shall  discuss  this  phenomenon  now. 
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3.3  Elimination  of  IPSPs  With  Maintenance  of 
Antagonistic  Center  Surround  Excitatory  In¬ 
put 

In  Section  2.4  we  discussed  ganglionic  models  in  which  the  interaction  be¬ 
tween  (non-DS)  excitatory  and  inhibitory  input  onto  the  ganglion  cell  un¬ 
derlie  the  cell’s  DS  response.  Also,  interpretation  of  ganglion  cell  membrane 
potentials  may  be  complicated  by  inhibitory  input  onto  the  cell  which  might 
not  have  a  direct  link  to  the  directional  properties  (Section  2.3.2). 

For  this  reason  we  have  recorded  from  ganglion  cells  using  whole  cell 
patch  electrodes  in  which  the  electrode  solution  is  free  of  ATP  (in  some 
cases  without  as  well)  (Appendix  B).  Given  the  large  bore  of  the 

electrode  (1-2  /im),  it  is  likely  that  the  cell  contents  are  dialyzed  by  the 
electrode  solution  within  several  minutes  after  the  st2irt  of  recording.  It  has 
been  reported  in  hippocampus  that  such  conditions  block  the  response  of 
GABAa  receptors  ([60]);  thus  this  technique  offers  a  method  for  selectively 
blocking  an  inhibitory  component  of  the  recorded  cell’s  s)rnaptic  input,  with¬ 
out  disturbing  the  test  of  the  network. 

Using  this  technique  we  have  recorded  from  ganglion  cells  which  show 
clear  light-evoked  IPSPs  and  EPSPs  at  normal  resting  potential  at  the  onset 
of  whole-cell  access.  However,  in  almost  all  these  cells,  within  typically  ten  to 
twenty  ml  .utes  after  whole-cell  act  err.,  I'-^ht-evoked  IPSPs  disappe  v;V  ' 
EPSPs  are  maintained,  suggesting  the  block  described  above.  This  apparent 
removal  of  direct  inhibitory  input  to  these  cells  wu  verified  by  depolarizing 
the  cells:  negative  synaptic  potentials  were  not  observed  despite  large  de¬ 
polarizations  (30  mV  above  the  resting  potential).  Hyperpolarizing  phases 
of  spike  repolarizations  were  observed,  however,  suggesting  preservation  of 
voltage-dependent  hyperpolarizing  mechanisms. 

Figures  B.6  through  B.12  in  Appendix  B  demonstrate  the  reduction  of 
IPSPs  to  flashing  spots  over  time.  In  this  chapter  we  shall  show  DS  record¬ 
ings  of  these  same  cells  or  cells  recorded  under  the  same  conditions.  To  re¬ 
capitulate,  DS  ganglion  cell  responses  without  measurable  inhibitory  input, 
especially  with  P/N  distinctions  that  are  stable  with  respect  to  hyperpolar- 
izing  or  depolarizing  current,  suggest  that  some  turtle  ganglion  cells  receive 
excitatory  input  which  is  already  DS. 
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3.4  Directionally  Selective  Responses 

We  shall  now  present  DS  recordings  from  several  gcmglion  cells  taken  with 
a  variety  of  protocols.  We  shall  describe  the  basic  findings  under  current 
clamp  and  voltage  clamp,  and  then  present  example  recordings  in  Figure  3.3 
through  3.25.  These  recordings  are  listed  in  Table  3.1  for  quick  reference. 

DS  data  is  presented  in  several  formats.  Polar  plots  of  single-valued 
response  parameters  are  presented  with  standard  deviations  when  the  pa¬ 
rameter  was  computed  for  each  trial,  and  then  averaged  over  several  trials. 
Response  waveforms  are  presented  both  in  the  standard  fashion  and  in  per¬ 
spective  plots  (moving  gratings  only).  Each  perspective  plot  is  composed  of 
a  series  of  averaged  current  or  voltage  clamp  waveforms  in  response  to  mo¬ 
tion  stimuli  at  16  orientations  (WF  protocol).  The  responses  shown  in  these 
plots  are  passed  through  a  low  pass  filter  (/e  =  20hz)  before  averaging.  As 
mentioned  earlier  in  the  description  of  moving  grating  burst  averaging,  each 
averaged  waveform  shows  both  the  transient  response  and  the  (beginning  of 
the)  steady  state  (grating)  response.  The  waveform  data  is  replicated  once 
by  concatenating  2  full  cycles  of  the  360®  sweep  of  the  stimuli  direction,  in 
order  to  more  clearly  show  any  directional  component. 

For  some  cells,  we  also  show  spot-mapped  receptive  fields  overlaid  with 
the  cells’  DS  axis. 

3.4.1  Directionally  Selective  Responses  Under  Current  Clamp 

DS  EPSPs  and  spike  responses  have  been  obtained  with  the  whole-cell  patch 
technique  under  the  conditions  which  reduce  or  eliminate  IPSPs.  The  DS  of 
the  responses  are  stable  with  respect  to  holding  current,  which  was  adjusted 
to  move  the  resting  potential  within  ±  30  mv  from  the  control  Erttt' 

3.4.2  Directionally  Selective  Responses  Under  Voltage  Clamp 

DS  responses  were  also  recorded  under  voltage  clamp.  In  the  perspective 
plots  of  these  responses  clamp  current  polarity  has  been  reversed  for  clarity; 
thus  positive  current  implies  an  excitatory  synaptic  event.  The  DS  of  the 
response  are  stable  with  respect  to  holding  potential,  which  ranged  from  -90 
to  -35  mV. 

For  example,  in  cell  A56  (Figure  3.8)  the  dependence  is  striking  (best  at 
about  157®),  and  is  independent  of  the  holding  potential  (Figure  3.10).  For 
comparison,  cell  A81  (Figure  3.13)  shows  little  dependence  of  the  voltage 
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clamp  waveforms  on  stimulus  direction. 

Although  the  degree  in  which  the  dendritic  tree  of  these  cells  is  actuadly 
clamped  during  the  voltage  clamp  protocol  is  not  know  directly,  the  obser¬ 
vation  of  DS  inward  currents  at  low  (somatic)  holding  potentials  is  further 
evidence  for  at  least  a  component  of  the  excitatory  synaptic  input  being  DS, 
independent  of  any  local  blockage  of  IPSPs. 
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CeU 

Type 

Figure 

Recording 

Method 

Stimulus 

Protocol 

Plot  Type 

A47 

± 

3.3 

CC 

WF3 

Conventional 

3.4 

CC 

WF3 

Polar 

VC 

WF15 

Conventional 

CC 

WF15 

Polar 

A56 

+ 

EC 

WF2 

Polar 

3.8 

VC 

WF4 

Perspective 

mtsM 

VC 

WF4 

Polar 

3.10 

VC 

WF6 

Perspective 

3.11 

VC 

WF6 

Polar 

3.12 

CC 

RFl-2 

RF  w/DS  axis 

A81 

± 

3.13 

VC 

WF6 

Non-DS  Persp. 

A84 

± 

3.14 

EC 

WF2 

Polar 

3.15 

CC 

WF3 

Conventional 

3.16 

VC 

WF5 

Perspective 

3.17 

VC 

WF5 

Polar 

A105 

4- 

3.19 

CC 

WF3 

Polar 

3.20 

CC 

PN2-4 

Conventional 

3.21 

CC 

RF2 

RF  w/DS  axis 

Alls 

+ 

3.22 

EC 

WF4 

Polar 

3.23 

CC 

PNl-3 

Conventional 

3.24 

CC 

WF5 

Polar 

3.25 

CC 

PN4-8 

Conventional 

3.26 

CC 

RF2 

RF  w/DS  axis 

Table  3.1:  Reference  to  figures  of  DS  recordings.  Cell  type:  ±  =  ON/OFF, 
+  =  ON.  Recording  Method:  CC  =  current  clamp,  VC  =  voltage  clamp,  EC 
=  extracellular.  Stimulus  Protocol:  WF  =  wide  field  DS  test  and  tuning,  RF 
=  receptive  field  map  with  spots,  PN  =  preferred/null  motion  runs.  Numbers 
after  stimulus  protocol  designations  refer  to  specific  experiments. 
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Figure  3.3:  DS  responses  to  moving  gratings  under  current  clamp  for  ON/OFF  ganglion 
cell  A47  (protocol  A47WF3).  Each  waveform  is  the  average  of  4  trials.  The  grating 
had  a  200  ^m  square  aperture,  a  temporal  frequency  of  1  hz  and  a  spatial  period  of 
800/zm.  In  this  protocol,  Op  =  67.5®,  i7/«um  =  0.2  (p  <  0.05)  (integral  of  V  - 
Er,it)-  This  recording  was  made  17  minutes  after  whole-cell  access  (there  was  no 
extracellular  response  recorded).  Although  this  cell  responded  to  both  ON  and  OFF 
edges  of  stationary  spots,  the  expected  2nd  harmonic  of  the  motion  response  is  small. 
Small  IPSPs  were  observed  within  5  minutes  of  whole-cell  access  during  receptive  field 
stationary  spot  mapping. 
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Figure  3.4:  Polar  plot  of  average  (V  -Ertit)iSD  (mV)  for  waveforms  whose  averages 
and  DS  statistics  are  shown  in  Figure  3.3  (A47WF3). 


Figure  3.5:  DS  responses  to  moving  grating  under  voltage  clamp  for  cell  A47  (protocol 
A47WF15),  taken  2  hours  after  the  responses  in  Figure  3.3.  Holding  potential  =  -80mV, 
average  of  4  trials.  Stimuli  are  identical  to  those  used  for  Figure  3.3.  Using  the  integral  of 
{IhoUing  —  Idamp)  **  the  response  variable,  the  DS  statistics  are  0p  =  112.5°,  Dlgum  = 
0.3  (p  <  0.05).  The  preferred  direction  has  apparently  shifted  from  the  previous  protocol 
(though  the  statistical  significance  of  this  was  not  evaluated),  but  remains  consistent 
with  the  original  axis  of  asymmetry. 


Figure  3.6:  Polar  plot  of  average  {IhoiSng  -  Idamp)  ^  SD  {pA)  for  waveforms  whose 
averages  are  shown  in  in  Figure  3.5  (A47WF15). 
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Figure  3.7:  Polar  plot  of  DS  extracellular  spike  response  of  ON  cell  A56  to  a  moving 
grating  (average  of  8  trials  per  orientation,  16  orientations).  Grating  aperture  was 
200;im  square,  with  a  temporal  frequency  of  1  hz  and  spatial  period  of  800/im.  For  this 
protocol  $p  =  135°,  DItytrn  =  1*0  (p  <  0.05).  The  standard  deviation  for  this  protocol 
was  not  available. 
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Figure  3.8:  Perspective  plot  of  DS  responses  to  moving  grating  under  voltage  clamp  for 
cell  A56,  taken  14  minutes  after  whole-cell  access  (average  of  2  trials).  The  response  to 
the  initial  ON  edge  has  a  delay  of  about  400ms,  and  a  very  slight  response  to  the  first 
OFF  edge  is  seen  around  1300ms.  The  holding  potential  =  -80mV  strongly  attenuates 
any  inhibitory  synaptic  current,  even  if  the  dendritic  tree  is  not  fully  clamped.  Thus, 
this  response  suggests  that  excitatory  input  to  this  cell  has  a  large  DS  component. 
Grating  aperture  was  lOO^m  square,  with  temporal  frequency  of  1  hz  and  spatial  period 
of  800/im. 
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Figure  3.9:  Polar  plot  of  average  {IhoUing  —  Idamp)  ±  SD  [pA)  for  waveforms  whose 
averages  are  shown  in  in  Figure  3.8  (A56WF4).  For  this  protocol,  9p  =  157.5°,  but 
there  is  no  statistically  significant  OS  based  on  this  parameter  alone  (presumably  because 
of  the  low  number  of  trials  ). 
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Figure  3.10:  Perspective  plot  of  DS  responses  to  moving  grating  under  voltage  clamp 
for  cell  A56,  taken  48  minutes  after  whoie-cell  access  (average  of  2  trials).  Stimulus 
parameters  identical  to  those  used  for  Figure  3.8.  Holding  potential  =  -35mV.  The 
lack  of  a  negative  phase  (relative  to  the  holding  current)  in  this  plot  of  the  inverted 
clamp  current  suggests  that  any  inhibition  activated  by  the  gratings  occur  in  phase  with 
the  excitation,  independent  of  orientation.  This  observation,  along  with  the  the  fact 
that  the  total  clamp  current  is  DS  suggests  that  either  excitation  or  inhibition  or  both 
types  of  synaptic  input  is  already  OS. 
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Figure  3.11:  Polar  plot  of  average  {Ihuding  —  Idamp)  ±  SD  (pA)  for  waveforms  whose 
averages  are  shown  in  in  Figure  3.10  (A56WF6).  As  for  the  A56WF4  data  shown 
in  Figure  3.8,  6p  =  157.5°,  but  there  is  no  statistically  significant  DS  based  on  this 
parameter  alone  (presumably  because  of  the  low  number  of  trials). 
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Figure  3.12:  Spot  mapped  receptive  field  for  cell  A56  (protocols  A56RF1  and  A56RF2). 
Grid  is  1000  by  1000  /im  (centered  on  electrode  tip),  spot  size  200  fjim  square,  flashing 
at  1  hz.  Spike  results  obtained  extracellularly;  EPSP  results  obtained  within  5  minutes 
of  whole-cell  access.  OS  orientation  from  A56WF  protocols  presented  in  earlier  fig¬ 
ures.  Circle  diameters  are  proportional  to  normalized  responses.  Both  spike  and  EPSP 
response  was  maximum  at  the  center  of  the  grid.  Note  that  the  OS  orientation  is  or¬ 
thogonal  to  the  somewhat  elongated  RF,  similar  to  results  reported  by  [28].  Also,  the 
EPSP  receptive  field  is  larger  than  the  spike  receptive  field,  as  would  be  expected  from 
a  thresholding  spike  mechanism. 
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Figure  3.13:  Perspective  plot  of  non-OS  responses  to  moving  grating  for  comparison 
with  DS  perspective  plots.  Voltage  clamp  protocol  for  cell  A81  (average  of  4  trials), 
holding  potential  =  — 80mV.  Grating  aperture  was  400  fim  square,  temporal  frequency 
of  2  hz.  This  cell  shows  no  dependence  on  the  stimulus  direction. 
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Figure  3.14:  Polar  plot  of  DS  extracellular  spike  responses  of  ON/OFF  ganglion  cell 
A84  to  a  moving  grating  (average  of  2  trials  per  orientation,  16  orientations).  Grating 
aperture  was  300/im  square,  temporal  frequency  of  2  hz.  DS  statistics  of  average  spike 
response  was  $p  =  112.5®,  Dl^rn  =  0.2  {p  <  0.05). 
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Figure  3.15:  DS  responses  to  moving  gratings  under  current  clamp  for  ON/OFF  gan¬ 
glion  ceil  A84  (protocol  A84WF3).  Each  waveform  is  the  average  of  4  trials.  The  grating 
had  a  200  /im  square  aperture,  a  temporal  frequency  of  2  hz  and  a  spatial  period  of 
800/im.  In  this  protocol,  dp  =  45®,  Dl^m  =  0.7  (p  <  0.05)  (spike  count).  This 
recording  was  made  9  minutes  after  whole-cell  access.  Note  that  the  ON  response  is 
much  weaker  than  the  OFF  response  under  current  clamp,  but  that  the  ON  EPSP  is 
still  slightly  directional. 


53 


1980  ms 


Figure  3.16:  Perspective  plot  of  response  to  moving  grating  under  voltage  clamp  for  cell 
A84,  taken  30  minutes  after  whole-cell  access.  Holding  potential  =  -80mV  (average  of 
4  trials),  suggesting  that  excitatory  input  to  this  cell  has  a  large  DS  component.  Grating 
aperture  was  200/im  square,  temporal  frequency  of  2  hz.  The  2nd  harmonic  from  the 
ON/OFF  cell  is  clearly  seen  in  this  motion  response,  and  the  ON  response  is  relatively 
larger  than  that  seen  in  Figure  3.15,  possibly  because  the  voltage  clamp  reduces  the 
contribution  of  inhibition,  and  the  inhibition  may  be  more  blocked  by  the  lack  of  ATP 
at  this  later  recording  time.  Also,  the  ON  response  adapts  more  than  the  OFF  response, 
which  in  turn  has  a  stronger  directional  dependence  over  time.  The  shift  in  preferred 
direction  between  this  figure  and  Figure  3.14  (here,  Bp  =  22.5°,  Dl^m  =  0.12  (p  < 
0.05))  might  be  accounted  for  by  a  shift  in  the  aperture  location  (see  Figure  3.18), 
although  the  statistics  of  a  broad  directional  tuning  may  also  account  for  an  apparent 
shift. 


Figure  3.17:  Polar  plot  of  average  {lading  -  hiamp)  ±  SD  {pA)  for  waveforms  whose 
averages  are  shown  in  in  Figure  3.16  (A84WF5). 


Figure  3.18:  Extracellular  determination  of  DS  properties  for  cell  A84  were  done  with 
reference  to  an  initial  estimate  (by  eye)  of  the  receptive  field  center  of  (z  =300, y  =0), 
units  in  ^m,  relative  to  the  electrode  tip.  Subsequent  and  more  careful  intracellular 
receptive  field  measurements  moved  the  receptive  field  center  to  (z  =350, y  =50). 
A  contribution  of  asymmetric  activation  of  the  symmetric  inhibitory  surround  in  the 
extracellular  OS  measurements,  with  the  DS  properties  of  correctly  centered  stimuli  can 
account  for  the  shift  in  $p  between  protocols  for  this  cell.  However,  since  the  extracellular 
9p  is  close  to  that  predicted  by  the  center/surround  DS  artifact,  it  is  possible  that  this 
DS  is  purely  due  to  center/surround  mechanisms. 
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Figure  3.19:  Polar  plot  of  DS  intracelluiar  spike  response  (with  SD)  of  ON/OFF  cell 
A105  to  a  moving  spot  taken  after  5  minutes  of  whole>cell  access.  Spot  is  200  /im  square, 
moving  at  4  ^m/msec,  with  a  path  length  of  1200/4m.  Average  of  2  trials  per  orientation, 
16  orientations.  The  OS  statistics  on  this  protocol  give  9p  =  45°,  Dl^im  =  0.47  (p  < 
0.05).  Although  the  strongest  responses  are  along  the  112.5/292.5°  axis  (showing 
a  strong  orientation  selectivity),  the  most  statistically  significant  radial  asymmetry  is 
towards  45  degrees. 
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Figure  3.20:  Single  trial  DS  responses  of  cell  A105  to  moving  spots,  7  minutes  after 
whole-cell  access,  with  positive,  0,  and  negative  holding  currents  (protocols  A105PN2- 
4).  Spot  stimulus  is  same  as  Figure  3.19,  except  that  orientation  is  restricted  along 
the  45/225”  axis.  Spot  motion  is  timed  so  that  the  spot  passes  over  the  center  of 
the  receptive  field  at  about  1000ms.  In  this  cell  the  preferred  orientation  of  45°  is 
independent  of  the  holding  current.  Also,  there  is  little  evidence  for  IPSPs  in  these 
traces,  suggesting  that  the  directionality  does  not  depend  on  ganglion  cell  inhibition. 
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Figure  3.21:  Spot  mapped  receptive  field  for  cell  A105  (protocol  A105RF2).  Grid  is  500 
by  500  nm  (center  at  (*100  /im,  300  fim)  relative  to  electrode  tip),  spot  size  200  ftm 
square,  flashing  at  1  hz,  responses  recorded  17  minutes  after  whole-cell  access.  OS 
orientation  from  A105WF  protocols  presented  in  earlier  figures. 
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Figure  3.22:  Polar  plot  of  DS  extracellular  spike  response  (with  SD)  of  ON  cell  A115 
to  a  moving  spot.  Spot  is  200/sm  square,  moving  at  2  ^m/msec,  with  a  path  length 
of  lOOO/im.  Average  of  2  trials  per  orientation,  16  orientations.  DS  statistics  give 
dp  =  270",  =  0.15  {p  <  0.05). 


Figure  3.23:  Single  trial  DS  responses  of  cell  A115  to  moving  spots,  immediately  after 
whole-cell  access,  with  positive,  0,  and  negative  holding  currents  (protocols  A115PN1- 
3).  Spot  stimulus  is  same  as  Figure  3.22,  except  that  orientation  is  restricted  along  the 
290/110°  axis  (timing  as  in  Figure  3.20).  Similar  to  cell  A105,  the  preferred  orientation 
of  290°  is  independent  of  the  holding  current,  and  there  is  little  evidence  for  IPSPs 
involvement  in  the  OS  response. 
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Figure  3.24:  Polar  plot  of  DS  intracellular  spike  response  (±517)  of  cell  A115  for  same 
protocol  as  Figure  3.22  (A115WF5).  This  protocol  was  taken  immediately  after  the 
P/N  runs  shown  in  Figure  3.23  (within  10  minutes  from  start  of  whole-cell  access).  The 
DS  statistics  for  this  protocol  give  9p  =  315®,  Dl^rn  =  0.22  (p  <  0.05). 
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Figure  3.25:  Responses  of  cell  A115  to  moving  gratings,  31  minutes  after  whole-cell 
access,  to  investigate  subunit  properties.  Each  waveform  is  the  average  of  2  trials, 
with  each  trial  consisting  of  2  full  cycles  of  the  grating.  Grating  apertures  are  200;xm 
square,  temporal  frequency  1  hz,  and  orientation  of  the  gratings  was  restricted  along  the 
290/110°  axis.  Aperture  centers  were  placed  at  100/im  steps  along  the  290/110°  axis, 
as  noted  in  the  figure.  In  this  series  local  DS  center/surround  responses  are  seen,  but 
the  DS  response  to  long  trajectory  stimuli  (such  as  in  Figure  3.23)  cannot  be  predicted 
from  the  localized  responses,  e.g.  there  is  no  consistent  DS  subunit  response. 
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Figure  3.26:  Spot  mapped  receptive  field  for  cell  A115  (protocol  A115RF1).  Grid  is  500 
by  500  fim  (center  at  (100  /xm,  -100  ftm)  relative  to  electrode  tip),  spot  size  100  /xm 
square,  flashing  at  1  hz,  extracellular  spike  response.  DS  orientation  from  A115WF 
protocols  presented  in  earlier  figures. 
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Chapter  4 

A  Model  of  Retinal 
Directional  Selectivity 


Models  for  retinal  directional  selectivity  have  proposed  various  locations  for 
the  asymmetric  path,  delay  elements,  and  non-linear  interactions  required 
for  this  computation,  as  we  discussed  in  Chapter  2.  The  experimental  results 
of  the  previous  chapter  add  some  constraints,  namely  that  all  these  elements 
must  occur  prior  to  the  DS  ganglion  cell,  which,  in  turn,  apparently  receives 
DS  excitatory  input. 

We  now  present  a  pre-ganglioni  -uodel  for  DS  whicK  ’ '  morphoiretricly 
similar  to  Vaney’s  ([66],  [64])  in  that  the  lateral  pathway  is  via  individual 
branches  of  amacrine  cells  with  tip  outputs  [12].  We  find  as  well  that  given  a 
plausible  set  of  constraints  on  the  distribution  and  biophysics  of  synaptic 
input  and  output  on  the  branches  and  b)  the  intrinsic  cable  properties,  the 
outputs  of  this  pathway  are  at  least  DD  and  normally  are  DS.  Thus,  both  the 
lateral  path  and  the  non-linear  interaction  underlying  the  DS  computation 
occur  on  the  same  substrate. 


4.1  Oriented  Amacrine  Cell  Branches  with  Di¬ 
rectionally  Selective  Outputs 

The  crucial  DS  element  in  this  circuit  is  shown  in  Figure  4.1,  and  the 
anatomy  of  the  asymmetric  pathways  and  the  target  DS  ganglion  cell  is 
shown  in  Figures  4.2  and  4.3.  We  suggest,  ta  Vaney,  that  DS  subunits  arise 
from  multiple  oriented  amacrine  cell  dendritic  tips  connecting  to  the  DS 
ganglion  cell.  In  the  context  of  retinal  DS,  the  directional  properties  of  a 
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cable  with  distributed  synaptic  conductance  input  has  also  been  described 
by  Grzywacz  and  Amthor  ([29]).  Likewise,  Koch,  Poggio  and  Torre  ([40], 
Figure  9a)  suggested  a  similar  (ganglionic)  arrangement  of  locally  symmet¬ 
ric  excitatory  and  delayed  inhibitory  inputs  along  the  branch  of  a  putative 
asymmetric  DS  ganglion  cell. 
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Receptive  Fields 

V  V  V 
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Figure  4.1:  (A)  Structure  of  cable  model  discussed  here.  In  this  model  the  time  element 
necessary  for  DS  is  provided  by  the  inherent  kine^'  of  the  input  synapr';*  (most  impor¬ 
tantly  by  the  slower  time  course  of  the  inhibition^,  and  the  nonlinearity  also  necessary 
for  OS  by  the  distributed  inhibitory  synaptic  conductance.  Under  certain  conditions  the 
location-dependent  delay  (phase  shift)  provided  by  the  cable  membrane  time  constant 
(Tm)i  and  an  output  nonlinearity,  may  also  be  relevant  (Sections  4.2.4  and  4.3.2).  The 
preferred  response  of  the  tip  output  above  is  for  motion  to  the  right,  for  the  model  ver¬ 
sion  analyzed  in  this  chapter.  (B)  Equivalent  circuit  for  a  section  of  the  dendrite  cable, 
including  the  axial  resistance,  the  linear  resting  membrane  resistance  and  capacitance, 
and  the  synaptic  inputs.  Directionality  of  the  distal  tip  output  arises  because  in  the 
null  direction  the  slower  inhibition  shunts  subsequent  excitatory  input  as  the  stimulus 
moves  away  from  the  tip.  In  the  preferred  direction  the  inhibition  is  much  less  effective 
at  attenuating  the  excitation,  due  to  the  on-the-path  condition  (see  text). 


66 


LIGHT 


Figure  4.2:  Relationship  between  the  asymmetric  pathways  in  the  cable  model  and 
the  target  DS  ganglion  cell,  and  pattern  of  activation  for  light  moving  in  the  preferred 
direction,  input  onto  the  asymmetric  amacrine  cell  dendritic  paths  are  from  symmetric 
(direct)  pathways  (either  via  bipolar  ceils  or  amacrine  cells).  The  DS  ganglion  cell  also 
receives  direct  input  from  bipolar  and  amacrine  cells  (not  shown).  Note  that  the  wave 
of  inhibition  lags  behind  the  wave  of  excitation.  Other  tips  of  the  amacrine  cell  with 
different  orientations  may  connect  to  other  DS  ganglion  cells,  since  the  tip  responses 
are  essentially  independent  of  one  another. 
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Figure  4.3:  Pattern  of  activation  for  light  moving  in  the  null  direction.  In  this  case, 
after  the  light  has  passed  over  the  amacrine  cell  dendrite  tips  there  is  a  zone  of  activated 
inhibition  between  the  tips  and  the  activated  excitation. 
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4.1.1  Start  by  Finding  the  Asymmetry 

Histological  data  from  retina  does  not  immediately  suggest  the  spatial  asym¬ 
metry  required  for  directional  selectivity,  as  mentioned  in  Chapter  2.  How¬ 
ever,  if  we  consider  putative  outputs  on  the  dendritic  tips  of  amacrine  cells, 
the  IPL  becomes  full  of  asymmetries,  since  for  a  large  number  of  amacrine 
cells  the  tips  are  displaced  with  respect  to  the  entire  tree  (see  Section  5.7). 
This  interpretation  deviates  from  the  assumption  of  classical  neuronal  pol2ur- 
ity:  input  onto  the  dendritic  tree,  conducted  proximally  and  integrated  by 
the  soma  membrane,  cell  output  via  the  somatic  axon.  Some  retinal  neurons 
challenge  this  dogma,  and  the  distinction  of  dendrite  versus  axon  becomes 
blurred.  Mciny  retinal  cells  do  not  have  an  axon  at  all  per  se,  and  single  cell 
processes  may  have  both  input  and  output. 

As  Vaney  pointed  out,  detailed  studies  of  starburst  amacrine  cell  den¬ 
drites  support  this  potential  asymmetry  for  the  rabbit  retina  ([26],  [25],  [65], 
[62]).  Small  swellings  in  the  outer  third  of  the  dendritic  tree  have  been  de¬ 
scribed,  and  electron  micrographs  show  that  synapses  with  ganglion  cells 
occur  in  the  swellings’  zone  [14].  In  contrast,  the  general  distribution  of 
S3rnaptic  inputs  along  the  dendrites  is  apparently  iiniform^. 

It  should  be  emphasized  that  the  critical  asymmetry  in  this  model  re¬ 
quires  inputs  all  over  the  amacrine  cell  dendrite  branch  and  outputs  near  the 
tip.  There  is  no  restriction  on  outputs  elsewhere  on  the  bramch,  however, 
ar. '  ae  most  gene~e.l  version  of  the  model  would  have  uniform  inputs  and 
outputs  over  the  entire  amacrine  cell  dendritic  tree.  On  the  other  hand,  only 
a  subset  of  those  outputs  will  have  strong  DS  properties,  and  an  assumption 
of  the  model  is  that  a  given  target  DS  ganglion  cell  will  receive  at  least 
some  of  the  outputs  from  this  subset.  This  problem  will  be  addressed  in 
Section  5.5. 

4.1.2  Add  Location  Dependence 

The  morphometric  asymmetry  is  necessary  for  directional  selectivity,  but 
not  sufficient.  If  the  dendrite  branch  is  isopotential  then  there  is  no  location 

’It  should  be  pointed  out  that  there  it  conflicting  support  for  the  starburst  amacrine 
cell  speciflcally  in  the  DS  circuit  for  rabbit.  For  example  Linn  and  Massey  [43]  suggest 
that  cholinergic  release  from  starburst  cells  is  not  inhibited  directly  by  GABA,  which  is 
contrary  to  the  model  prediction.  However,  other  amacrine  cells  may  have  the  necessary 
input/output  relationship.  Later  we  shall  show  simulations  of  a  cell  whose  morphometry 
is  based  on  the  starburst  cell;  despite  the  evidence  against  the  starburst  ceU  speciflcally, 
we  feel  that  the  conclusions  from  this  tree  geometry  are  applicable  to  other  amacrine  cells. 
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dependence  of  branch  input,  precluding  any  sort  of  motion  selectivity,  much 
less  directional  selectivity.  When  we  add  the  intracellular  resistance  of  the 
cable,  however,  then  location  dependence  results. 

4.1.3  Directional  Dependence  for  the  Linear  Case 

Let  us  now  examine  the  response  of  the  cable  tip  to  moving  inputs.  We 
begin  with  the  simplest  case:  distributed  excitatory  sjmaptic  input  onto 
the  (linear)  cable.  If  the  synapses  are  modelled  as  current  sources,  e.g. 
the  excitatory  synapse  injects  a  depolarizing  current,  then  the  tip  response 
can  be  obtained  by  the  appropriate  superposition  of  responses  to  individual 
inputs.  The  resulting  waveforms  for  motion  in  the  two  directions  will  have 
different  shapes,  but,  since  the  system  is  linear,  equal  areas. 

As  simulations  by  Rail  demonstrated  implicitly  ('^6]),  we  now  have  in 
the  dendrite  cable  the  prerequisites  -  location  dependen  :y  and  asymmetry  - 
for  DD  (but  not  DS).  Accentuating  the  difference  is  the  location-dependent 
delay  of  the  cable  (due  to  the  cable  time  constant).  This  delay  tends  to  in¬ 
crease  the  output  amplitude  for  motion  towards  the  tip;  the  EPSPs  arriving 
at  the  tip  are  more  correlated  since  the  cable  delay  (phase  shift)  counters 
the  motion  delay.  Otherwise  the  two  waveforms  would  simply  be  reversed 
in  time.  We  shall  return  to  this  phenomena  later  in  the  simiilations  section 
(Section  4.3.2). 

Now  •<  -me  that  at  each  point  on  the  cable  there  is  a  paired  eici- 
tatory/inhibitory  input,  with  the  inhibitory  kinetics  slower  than  that  of  the 
excitation^.  Again,  if  the  inhibitory  synapse  is  modelled  as  a  hyperpolar- 
izing  current  sotirce  the  (linear)  responses  for  the  two  directions  will  have 
equal  area.  The  DD  response  holds,  still  without  DS. 

4.1.4  Synaptic  Nonlinearities  and  Cable  Directionality 

Synaptic  input,  however,  is  more  accurately  modelled  by  a  conductance 
change  in  series  with  a  battery.  The  circuit  is  now  a  nonlinear  one:  the 
variable  conductance  of  the  inputs  means  that  motion  responses  cannot  be 
derived  from  the  superposition  of  a  sequence  of  point  responses.  As  a  result 
the  integral  of  the  response  for  the  two  directions  is  generally  not  conserved. 

The  effect  of  the  inhibitory  input  is  particularly  important  on  the  direc¬ 
tionality  of  the  cable  tip  response.  Assuming  that  the  reversal  potential  for 

^Note  that  if  the  excitatory  and  inhibitory  kiaetict  were  identical,  then  eynapie*  at  the 
tame  location  could  be  lumped  together  at  one. 
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the  inhibition  is  close  to  the  cell's  resting  potential,  then  the  nuun  effect  of 
the  inhibition  will  be  to  shunt  locally  any  excitatory  current  to  the  extracel¬ 
lular  space.  The  effect  of  the  shunt  on  the  tip  response  strongly  depends  on 
the  relative  location  of  the  excitatory  and  inhibitory  inputs  with  respect  to 
the  output.  As  Hall  ([56])  and  Koch  et  al.  ([41])  showed,  inhibition  is  most 
effective  when  it  is  on  the  path^  that  is  interposed  between  excitation  and 
the  output.  Conversely,  if  the  excitation  is  closer  to  the  output,  then  the 
inhibitory  shunting  is  much  less  effective  (Figure  4.4). 

The  inhibitory  shunting  now  supplies  the  necessary  nonlinearity  for  the 
DS  response.  This  prediction  is  in  concert  with  the  experimental  evidence  as 
to  the  importance  of  inhibitory  mechanisms  for  DS,  and  with  the  evidence 
that  inhibitory  shunting  might  mediate  DS  in  rabhit  ([3]).  We  also  note  that 
while  the  efficacy  of  the  inhibitory  synapse  with  respect  to  DS  is  mainly  due 
to  its  shunting  component,  it  is  also  true  that  the  precise  value  for  the 
inhibition  Ertv  is  not  crucial  (as  long  as  the  Ervo  is  in  the  neighborhood 
of  Locality  of  synaptic  interaction  is  strongest  when  the  inhibitory 

Ertv  =  Ere$t  (silent  inhibition),  but  this  is  not  relevant  for  the  geometry  of 
the  model  presented  here. 

4.1.5  Considering  Tip  Output  Nonlinearity,  Facilitation,  and 
the  Sign  of  the  Output 

Inhibition  onto  thv.  cw'.—  suffices  to  make  tlie  towards-tip  response  larger 
than  the  opposite  direction,  so  a  linear  output  function  at  the  tip  would 
preserve  the  distinction.  It  is  more  likely,  though,  that  the  output  synapse 
has  a  threshold.  However,  the  resulting  supralinear  region  will  only  amplify 
the  DS  distinction.  Satiuation  of  the  output,  on  the  other  hand,  may  have 
other  consequences  (Section  4.2.4). 

At  first  glance,  the  experimental  evidence  of  preferred-direction  facilita¬ 
tion  data  mentioned  earlier  supports  placing  a  facilitatory  mechanism  on  the 
DS  pathway,  as  opposed  to  only  a  symmetric  location  (e.g.  intrinsic  to  the 
ganglion  cell  spike  mechanism).  Further  evidence  for  a  DS  path  location  for 
the  facilitatory  mechanism  in  rabbit  ON/OFF  DS  cells  has  been  reported  by 
Grzywacz  and  Amthor  ([31]),  including  ON  and  OFF  segregation  of  facilita¬ 
tion,  no  correlation  between  facilitation  strength  and  response  strength  as  a 
function  of  receptive  field  position,  and  no  correlation  between  time  courses 
of  f^ilitation  and  response.  These  results  could  be  accounted  for  by  vari¬ 
ous  mechanisms  located  at  the  aznacrine  cell  tip.  For  example,  facilitatory 
mechanisms  intrinsic  to  the  synaptic  output  might  suffice.  Likewise,  we  may 
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Figiire  4.4:  Th«  effect  of  inhi'oii'  ''n  or  tin  output  is  dependent  on  the  relative  positions 
of  the  excitation  and  inhibition  with  respect  to  the  output.  To  illustrate,  we  use  an 
extreme  case  in  which  the  inhibitory  battery  is  equal  to  the  resting  potential,  and  the 
activated  inhibitory  conductance  is  infinite.  In  (A)  the  excitatory  input  is  closer  to  the 
tip  than  the  inhibitory  input.  Although  the  membrane  potential  at  the  site  of  inhibition 
is  clamped  to  the  resting  potential  (short-circuited),  the  axial  resistance  of  the  cable 
partially  isolates  the  excitatory  location,  and  the  tip  (V^)  is  depolarized.  In  (B)  the 
interposed  inhibitory  input  clamps  the  cable  between  the  excitatory  input  and  the  output 
to  the  resting  potential,  so  that  the  output  stays  at  rest.  Adding  membrane  capacitance, 
a  finite  conductance  to  the  activated  inhibitory  synapse,  or  making  the  inhibitory  reversal 
potential  not  exactly  equal  to  the  resting  potential  does  not  change  the  basic  interaction. 
However,  a  non-zero  RaxUd  **  crucial. 
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consider  a  voltage  and  time-dependent  K'*'  channel  at  the  tip  which  is  nor¬ 
mally  open  at  rest  but  slowly  inactivates  with  depolarization,  similar  to  the 
Id  channel  identified  in  hippocampal  pyramidal  cells  (see  review  by  Storm 
[61]).  The  characteristics  of  this  channel  would  cause  the  tip  output  to  be 
“primed”  by  distant  cable  input,  such  that  subsequent  excitation  near  the 
tip  (for  motion  in  the  preferred  direction)  would  be  unopposed  by  the  now- 
inactivated  K'*'  shunt.  In  Section  4.3.3  we  show  simulations  of  an  apparent 
motion  protocol  that  demonstrates  this  particular  facilitatory  mechanism. 

Two  versions  of  this  model  include  either  that  the  amacrine  DS  output 
forms  an  excitatory  connection  with  a  DS  ganglion  cell,  or  the  connection  is 
inhibitory.  It  is  also  possible  to  include  plausible  voltage  and  time  dependent 
nonlinearities  in  the  cable  so  that  the  preferred  direction  for  the  branch 
output  is  away  from  the  tip.  Also,  the  tip  output  may  pass  through  sign- 
inverting  bipolar  cells  or  amacrine  cells  before  the  ganglion  cell. 

While  any  combination  of  the  above  polarities  will  yield  a  ganglion  cell 
DS  response,  the  fact  that  directional  responses  have  been  recorded  in  the 
absence  of  GABAergic  inhibition  (Section  2.3.1)  suggests  that  the  DS  path¬ 
way  has  an  excitatory  component  at  every  junction.  Further,  the  data  de¬ 
scribed  in  Chapter  3,  in  particular  with  local  block  of  GABAergic  inhibitory 
input,  supports  the  excitatory  connectivity  version.  Also,  simulations  of 
model  cells  with  likely  membrane  parameters  favor  the  preferred-towards- 
tip  orientation  for  the  cable  output  (Section  ^.3  and  unpublished  data). 


4.2  Predictions  of  the  Model 

We  now  consider  testable  predictions,  namely  that  for  somatic  recordings  of 
amacrine  and  target  ganglion  cells. 

4.2.1  Directionally  Selective  Somatic  Recordings? 

With  respect  to  amacrine  cells,  the  model  predicts  that  moving  stimuli  cen¬ 
tered  on  the  (somatic)  receptive  field  wotild  not  elicit  a  directionality,  as¬ 
suming  a  symmetrical  dendritic  tree.  If  either  the  stimuli  or  the  dendritic 
tree  was  asymmetric  with  respect  to  the  soma,  then  a  directional  response 
would  result,  perhaps  similar  to  the  preganghonic  DS/DD  recordings  of  De- 
Voe  et.  al.  ([23]).  While  this  result  supports  the  proposed  mechanism  for 
DS,  it  does  not  link  the  mechanism  to  an  identified  DS  ganglion  celL 

From  the  point  of  view  of  DS  ganglion  ceil,  the  modd  predicts  that  the 
input  waveforms  to  the  cell  are  themselves  DS;  any  ganglionic  computational 
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mechanisms  will  be  inherently  symmetric  and  serve  only  to  refine  the  P/N 
response  properties.  This  result  is  in  contrast  with  the  ganglionic  model 
class,  in  which  the  inputs  are  not  DS  by  themselves. 

4.2.2  Dependence  of  Directional  Selectivity  on  Ganglion  Cell 
Membrane  Potential 

Differences  in  P/N  EPSPs  are  predicted  by  the  ganglionic  model,  but  since 
the  differences  arise  from  ganglionic  interaction  of  conductances  with  dif¬ 
ferent  reversal  potentials,  this  model  predicts  that  the  relationship  between 
the  preferred  and  null  EPSPs  will  depend  on  the  membrane  potential. 

For  example,  let  us  consider  the  ganglionic  AND-NOT  circuit:  Null  di¬ 
rection  response  reflects  a  temporal  overlap  of  the  excitatory  input  with 
inhibitory  (mainly  shunting)  input  at  the  ganglion  cell.  Since  the  inhibitory 
reversal  potential  might  be  near  the  rest  potential,  there  would  be  little  neg¬ 
ative  portions  in  either  the  preferred  or  null  response.  A  negative  portion 
would  result  if  the  membrane  potential  is  artificially  raised  by  injecting  cur¬ 
rent.  The  “unmasked”  IPSP  would  be  expected  to  be  more  correlated  with 
the  control  EPSP  in  the  null  direction  as  compared  to  the  preferred  direction 
(e.g.  [47]).  If  the  ganglion  cell  was  hyperpolarised  by  the  electrode,  the  “un¬ 
masked”  inhibitory  input  would  now  contribute  a  component  to  the  EPSP. 
With  sufficient  hyper /depolarization  the  amplitude  of  the  null  response  could 
become  greater  than  that  of  the  preferred  response  (see  i  i..l  and  2.2). 

This  result  is  in  clear  distinction  to  a  pre-ganglionic  model.  The  single 
reversal  potential  of  the  ganglion  cell  DS  inputs  imply  that  the  ratio  of 
the  preferred  and  null  EPSPs  is  independent  of  the  membrane  potential: 
Manipulating  the  membrane  potential  will  not  reverse  the  P/N  axis. 

4.2.3  Comparing  Total  Synaptic  Input  for  the  Preferred /Null 
Responses 

The  above  result  is  related  to  measuring  the  somatic  input  conductance, 
Gin{S,  t),  of  the  ganglion  cell  during  a  motion  stimulus  5(e,  t).  For  a  lumped 
cell  approximation  and  no  voltage-dependent  membrane,  predictions  about 
GiniS,t)  are  simple.  In  the  ganglionic  model,  the  only  difference  in  the 
inputs  for  the  preferred  versus  null  responses  is  in  their  timing.  Thus  the 
model  predicts: 


J  GUS{x,t))  =  y  G/n(5(-«,f)) 
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where  S{x,t)  and  S{—x,t)  are  stimuli  moving  in  the  preferred  and  null 
directions,  respectively.  On  the  other  hand,  for  the  pre-g{inglionic  model 
there  is  more  synaptic  input  for  one  direction  versus  the  other;  therefore 
(asstiming  the  lumped  cell  without  voltage  dependencies): 

jGiniSix,t))>  I GUSi-x,t)) 

for  a  pre-gimglionic  excitatory  DS  input  model,  and 

J  Gtn(S(x,t))  <  I  Gln{S{-X,t)) 

for  a  pre-ganglionic  inhibitory  DS  input  model.  An  advantage  of  comparing 
the  integrals  of  (?;„  for  the  preferred  and  null  stimuli  is  that  the  effect  of 
the  membrane  capacitance  may  be  ignored.  It  can  be  shown  that  under 
some  constraints  on  the  ganglion  ceil  and  the  experimental  protocol,  similar 
relationships  are  testable  even  with  distributed  inputs  on  dendritic  trees  and 
voltage  dependent  membrane. 

4.2.4  Dynamic  Range  of  Cable  Mechanisms:  Saturation  and 
Reversal  of  Directional  Selectivity 

Real  biophysical  mechanisms  saturate,  e.g.  a  supralinear  transfer  function 
will  not  stry  supralineer  for  unbounded  inputs.  We  now  consider  possible 
implicatioxib  of  saturation  on  the  performance  of  the  model. 

Saturation  of  an  input  synaptic  conductance  onto  the  cable  will  not 
change  the  basic  distinctions  between  the  preferred  and  null  waveforms,  since 
this  saturation  is  independent  of  stimulus  direction.  Saturation  (strictly 
speaking,  a  sublinear  region)  of  the  output  nonlinearity  can  have  quite  dif¬ 
ferent  effects:  for  strong  enough  cable  excitation  the  null  waveform,  with 
its  greater  temporal  support  (due  to  cable  Tn  spreading  out  the  synaptic 
potentials),  will  eventually  yield  a  final  output  that  is  equal  to  the  preferred 
output,  despite  the  larger  amplitude  of  the  preferred  waveform.  Increasing 
excitation  further  could  cause  a  reversal  of  the  P/N  orientation,  after  the 
output  nonlinearity.  This  might  be  observed,  for  example,  with  a  stimulus 
contrast  that  is  normally  DS,  but  under  conditions  in  which  inhibition  is 
reduced  or  blocked  (e.g.  with  pharmacological  manipulations).  Whether 
DS  would  be  entirely  eliminated  or  reversed  would  be  dependent  on  circuit 
and  stimulus  parameters. 

This  result  is  unique  to  models  in  which  the  DS  asymmetry  includes 
distributed  excitatory  input  with  relative  delays  along  the  P/N  axis.  Such 
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models  been  explored  theoretically  for  the  interpretation  of  DS  reversal  in 
fly  ([54]),  amd  in  Section  4.3.2  we  show  simulations  which  demonstrate  this 
effect. 


4.3  Simulations  of  Morphometricly  and  Biophys- 
ically  Detailed  Amacrine  Cell  Models 

To  investigate  neuronal  properties  which  are  pertinent  to  the  cable  model, 
we  have  run  simulations  of  amacrine  cells.  The  model  parameters  are  as  con¬ 
strained  as  possible;  morphometry  is  obtained  from  histological  data  in  the 
literature,  and  membrane  properties  are  either  inferred  from  experimental 
data  and/or  supported  by  theoretical  studies  of  other  cells  ([9]). 

An  important  aspect  of  these  simulations  is  that  dynamic  retinotopic 
stimuli  may  be  used  as  input  to  the  model  circuit.  Since  the  simulator 
maintains  the  3-dimensional  structure  of  the  cells,  it  is  straightforward  to 
interpret  model  response  to  realistic  arbitrary  stimuli  (Appendix  C). 

4.3.1  Simulations  of  Asymmetric  Responses  from  Symmetric 
Cells 

In  Figures  4.5  through  4.7  we  show  a  simulation  of  an  amacrine  cell  whose 
morphometry  is  taken  ..m  a  rabbit  starburst  sunacrine  cell  ([62]).  A  U  mem¬ 
brane  and  cable  properties  are  uniform,  including  specific  capacitance  (Cm), 
membrane  resistivity  (Am))  intracellular  resistivity  {Ri),  excitatory  synaptic 
density  (Cec),  and  inhibitory  synaptic  density  (Cm).  Am  and  At  are  fixed 
at  lOGAflcm^  and  200(lcm  ([58]),  respectively.  Cm  is  set  to  1.0  /iF/cm^ 
(to,  =  lOOmS),  Cm  is  set  to  100pS//im^,  and  Ces  to  IpS/^m^.  The  resting 
potential  for  the  cell  is  -70mV.  The  reversal  potentials  for  the  excitatory 
and  inhibitory  synapses  are  50mV  and  -70mV,  respectively.  No  voltage- 
dependent  membrane  is  included. 

The  synaptic  transfer  functions  are  fixed.  The  excitatory  input  response 
is  given  by  the  half-wave  rectified  convolution  of  the  light  stimulus  by  the 
difference  of  two  alpha  functions,  Or(t), 


the  first  with  a  r  of  10  mS  and  the  second  with  a  r  of  60  mS.  This  transient 
response  is  t}rpical  of  retinal  ON  response  to  a  flashing  stimulus.  The  in¬ 
hibitory  response  is  the  half-wave  rectified  convolution  of  the  light  stimulus 
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with  a  single  alpha  function  of  unit  area  and  r  =  100  mS.  This  sustained 
ON  response  approximates  the  time  course  for  inhibition  to  various  stimuli 
[2].  Adding  am  OFF  component  to  the  transfer  functions  does  not  change 
the  motion-dependent  behavior.  For  simplicity  the  spatial  impulse  response 
for  both  synapse  types  is  although  the  receptive  field  of  candidate 

bipolar  or  amacrine  cell  outputs  are  on  the  order  of  tens  to  hundreds  of 
microns  wide.  Since  spatial  properties  derive  from  the  successive  convolu¬ 
tion  of  each  circuit  element,  minimizing  the  intrinsic  spatial  spread  of  the 
amacrine  input  yields  an  upper  bound  on  overall  spatial  resolution  of  the 
DS  dendrite  structiire. 

The  large  value  for  is  consistent  with  values  measured  in  our  lab  and 
others  in  various  neurons  using  the  whole  cell  patch  technique  (e.g.  [19]). 
Synaptic  conductance  densities  and  kinetics  are  less  constrained  as  far  as  the 
literature  is  concerned;  the  range  of  values  we  have  chosen  produce  synaptic 
potentials  which  are  consistent  with  available  data. 

In  Figure  4.6  we  show  the  response  of  the  soma  and  a  distal  node  to  a 
bar  moving  across  the  entire  field  of  the  ceU.  The  soma  response  is  almost 
identical  for  opposite  motions,  while  the  tip  response  is  highly  DS.  In  Figure 
4.7  we  have  plotted  the  total  response  of  each  dendritic  tip  of  the  cell,  scaled 
by  stimulus  speed,  as  a  function  of  angle  off  the  soma,  for  a  slit  traveling 
across  the  entire  breadth  of  the  cell.  Despite  the  overall  symmetry  of  the 
cell,  this  functional  is  highly  directional. 

4.3.2  Parametric  Simulations  of  Cable  Mechanisms 

In  this  section  we  present  a  series  of  simulations  on  a  simple  cell  model 
(chosen  to  allow  rapid  parametric  simulations)  to  illustrate  how  the  different 
mechanisms  in  the  dendrite  cable  interact.  The  basic  structure  of  this  model 
is  a  1-dimensional  symmetrical  cell  with  two  opposing  unbranched  processes 
(Figure  4.8).  Cell  membrane  and  S3maptic  parameters  are  identical  to  the 
simulated  cell  in  Figure  4.5,  except  as  follows:  We  shzdl  vary  Cm  (1.0  and 
1.0~*  /iF/cm*,  or  =  100ms  and  0.001  ms)®,  and  Gin  (100,  10,  1,  0.1, 
and  0  pS//im®).  The  goal  here  is  to  see  how  directionality  depends  on  the 
inhibitory  input  versus  the  intrinsic  cable  delay  properties. 

As  before,  the  stimxilus  is  a  moving  slit,  50^m  wide,  now  with  velocities 
of  0.5,  1,  2,  4,  8,  and  16  /im/msec.  The  response  waveforms  from  the 
righthand  distal  tip  for  rightward  and  leftward  motion  will  be  compared. 

*A  Cm  of  0  itF/cBO?  wu  not  possible  due  to  the  integration  technique  of  the  simulator. 
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Figure  4.5:  Simulations  were  done  of  a  model  cell  based  of  the  morphometry  of  a  rabbit 
starburst  amacrine  cell,  taken  from  Tauchi  and  Masland  ([52]).  Shown  here  is  a  flat 
mount  view  of  cell,  with  snapshots  of  the  stimulus  slit  trajectory.  Node  32211  referred 
to  in  Figure  4.6  is  the  distal  tip  farthest  to  the  right. 


In  particular,  we  shall  compute  a  directional  index  (DI)  (after  [32])  for 
both  the  linear  integral  of  the  waveforms  and  the  integral  of  the  waveforms 
after  being  passed  through  a  sigmoidal  nonlinefirity  (representing  synaptic 
transmission): 


p,  I  f(Vp)  -  S  fjVN) 

ff{Vp)  +  ff{VN) 


DI=  -I  £>/  =  0  DI=1 

Strongly  Selective  for  V/v  No  Directionality  Strongly  Selective  for  Vp 

(Motion  Away  From  Tip)  (Motion  Towards  Tip) 

where  /()  is  either  the  identity  (=>  DJ{Average))^  or  the  sigmoidal  nonlin¬ 
earity  shown  in  Figure  4.12  (^  DI(Sigmoid  Average)).  Vp  is  the  voltage 
waveform  for  the  distal  tip  for  light  moving  toward  the  tip,  and  Vtf  is  the 
waveform  for  the  opposite  direction. 

As  shown  in  Figtires  4.11  and  4.12,  the  “normal”  directionality  of  the  tip 
outputs  {DI  >  0)  is  strongly  dependent  on  the  presence  of  inhibition.  How¬ 
ever,  if  a  nonlinear  functional  is  applied  to  the  tip  waveform,  the  distributed 
nature  of  the  cable  without  inhibition  generates  a  directional  response,  al- 
thcugh  it  is  much  weaker  than  the  control  case.  In  addition,  the  particidar 
nonlinearity  and  other  model  parameters  in  Figure  4.12  cau.,  3s  a  DI  reversal 
when  inhibition  is  blocked.  As  illustrated  in  Figures  4.9  and  4.12,  the  ability 
of  a  nonlinear  integrator  to  nuike  this  distinction  depends  on  the  delay  pro¬ 
vided  by  the  cable  capacitance,  since  a  time  delay  is  a  necessary  condition 
for  DS  (or  DD).  This  prediction  is  consistent  with  the  experimental  results 
cited  earlier  for  turtle  and  fly. 

4.3.3  Simulations  of  Facilitation  by  an  Inactivating  K'^  Chan¬ 
nel 

Apparent  motion  protocols  consisting  of  two  stimuli  in  the  sequence  X  — »  T 
can  be  used  to  evaluate  how  the  relative  position  of  X  and  Y  affects  the 
response  to  Y,  compared  to  the  response  to  Y  alone.  In  a  DS  system,  either 
X  must  reduce  the  Y  response,  when  the  sequence  is  in  the  null  direction, 
and/or  X  must  enhance  the  Y  response  for  the  preferred  direction  sequence. 

As  mentioned  in  Chapter  2,  both  the  former  (null  direction  inhibition)  and 
the  latter  (preferred  direction  facilitation)  have  been  observed. 
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In  Figiires  4.13  and  4.14  we  show  a  simulation  of  apparent  motion  for 
the  starburst  amacrine  cell  presented  earlier.  The  cell  parameters  are  iden¬ 
tical  to  before,  except  that  an  inactivating  K'^  channel  (as  mentioned  in 
Section  4.1.5)  wa^  included  in  the  distal  portions  of  the  tree.  In  Figure  4.14 
the  response  of  the  second  stimulus  is  strongly  affected  for  motion  in  both 
directions. 

4.3.4  Accounting  for  the  Anomalous  Area  Response 

As  mentioned  in  Section  2.4  Barlow  and  Levick  reported  a  subregion  of  the 
DS  ganglion  cell's  spot  receptive  field  in  which  preferred  and  null  responses 
were  similar.  This  subregion  lay  on  the  spot  receptive  field  edge  firom  which 
preferred  motion  originated.  Since  both  directions  give  a  response,  the  sim¬ 
plest  explanation  is  that  the  receptive  field  of  the  DS  sub- circuit  also  has 
an  ’’anomalous  area”,  which  is  recapitulated  by  the  ganglion  cell.  For  a  rea¬ 
sonable  set  of  parameters  our  model  predicts  such  an  area  for  a  particular 
morphological  cletss  of  amacrine  cells:  those  amacrine  cells  whose  individual 
branches  have  orthogonal  mates  (e.g.  stellate  and  starburst  morphological 
types). 

Consider  the  response  of  a  tip  output  of  such  a  cell  to  the  following 
stimulus:  a  moving  bar  whose  length  is  greater  than  the  threshold  spot 
diameter  for  motion,  and  whose  trajectory  is  over  some  fraction  of  the  cell’s 

it'itic  field,  along  the  preferred/  null  axis  of  the  dp.  Movement  o/er 
the  distal  third  of  the  tree  adjacent  to  the  output  gives  the  expected  DS 
waveform.  Movement  over  the  distal  third  on  the  other  side  of  the  soma 
produces  a  much  weaker,  perhaps  nil,  response,  due  to  the  distance  and  the 
intervening  somatic  load.  However,  movement  centered  over  the  soma  would 
produce  a  more  nearly  equivocal  set  of  preferred  and  null  responses.  This 
is  because  the  bar  passes  activates  synapses  on  dendrites  with  a  trajectory 
that  is  orthogonal  to  those  dendrites’  preferred/nuU  axes.  These  branches 
get  activated  almost  all  at  once  as  the  bar  sweeps  over  the  soma.  The 
resulting  bolus  of  excitatory  current  might  by  strong  enough  to  generate  a 
tip  response,  despite  the  cable  attenuation. 

4.3.5  An  Asymmetric  Substrate  for  Non-Linear  Interactions 
Avoids  Preferred  Direction  Inhibition 

As  described  in  Section  2.4.1,  putting  a  S3maptic  vetoing  mechanism  onto  a 
symzaetric  ganglion  cell  dendritic  tree  has  the  problem  of  preferred  direction 
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inhibition  for  motion  over  the  side  of  the  receptive  field  that  lies  towards 
the  preferred  direction.  In  these  models  there  is  a  tradeoff  that  depends  on 
the  duration  of  the  inhibitory  response:  the  longer  the  response,  the  wider 
the  null  velocity  bandwidth,  but  also  the  greater  the  preferred  direction 
inhibition.  The  model  described  here  does  not  have  this  tradeoff,  since  the 
integrating  substrate  (the  single  amacrine  cell  dendrite)  is  asymmetric.  The 
important  signal  flow  is  always  in  one  direction  (towards  the  dendrite  tip). 
For  a  single  spot,  there  is  no  velocity  nor  duration  of  inhibitory  kinetics  such 
that  preferred  direction  motion  will  elicit  residual  (interposed)  inhibition 
that  can  interfere  with  excitation. 

4.4  Summary 

We  have  presented  a  model  for  directional  selectivity  in  retina  whose  anatom¬ 
ical  structure  and  simulated  performance  is  consistent  with  the  data.  The 
crucial  elements  for  directionality  are  as  follows: 

•  As3rmmetric  input  to  output  distribution  on  the  DS  amacrine  cell  den¬ 
drite  pathway. 

•  Intracellular  resistivity  of  the  dendrite  cable  allowing  on — the-path 
interactions  which,  in  turn,  depend  on  stimulus  direction. 

•  Inhibitory  nonlinear  shunting  of  sufficient  duration  to  mask  subsequent 
proximal  excitation  of  the  distal  tip  output. 

Further,  we  have  showed  that  other  cable  mechanisms  may  generate 
testable  predictions  unique  to  this  model. 
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Figtire  4.G:  Responses  of  cell  in  Figure  4.5  to  50/im  slit  moving  from  left  to  right  and 
back  again  at  2  ^m/ms.  The  voltage  at  the  rightmost  distal  tip  in  Figure  4.5  (Node 
32211)  is  highly  DS,  whereas  the  soma  response  is  barely  DD.  Model  parameters  are 
given  in  the  text.  Simulations  also  show  that  inclusion  of  a  conductance  which 
inactivates  with  depolarization  can  both  facilitate  the  preferred-direction  response  and 
attenuate  the  null-direction  response  (see  Sections  4.1.5  and  4.3.3). 


«  T^p  Int€flrt1  (riU»<«cond»)  Sc«1»d  by  Stinulup  Speed  (uH/wi 11 i«»cond)  I 

Figure  4.7:  Polar  plot  of  integral  of  distal  tip  voltages  scaled  by  the  stimulus  speed 
versus  angle  of  tip  with  respect  to  the  soma,  relative  to  the  slit  trajectory  shown  in 
Figure  4.5.  Slit  speeds  include  0.5,  2,  and  8  /zm/ms.  Despite  the  overall  symmetry  of 
the  cell,  the  tip  outputs  respond  asymmetrically  to  motion. 
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Figrue  4.8:  Model  cell  with  two  opposing  dendrites  whose  diameters  taper  linearly  from 
1.0  to  0.2  fim,  (proximal  to  distal),  with  snapshots  of  the  stimulus  slit  trajectory.  The 
output  used  in  the  later  figures  is  from  the  righthand  distal  tip. 
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Figure  4.9:  Example  set  of  simulations  for  cell  in  Figure  4.8  for  a  stimulus  moving  at 
4  ^im/msec  at  various  values  for  Gin  and  Cm-  On  the  left  Cm  =  1.0/xF/cm^,  and  on 
the  right  Cm  =  1<0~^  /xF/cm^.  Tip  output  is  from  the  righthand  distal  tip  in  Figure  4.8, 
and  the  solid  and  dashed  waveforms  are  for  leftward  and  rightward  motion,  respectively. 
When  inhibition  is  absent  (bottom  simulations)  the  cable  capacitance  distorts  the  equal 
area  waveforms  so  that  a  nonlinear  integrator  may  distinguish  them  (bottom  left).  If  the 
nonlinearity  saturates,  then  the  “fatter”  leftward  response  yields  a  larger  total  output 
than  the  higher  amplitude  rightward  response  (DS  reversal).  With  no  inhibition  and 
very  low  capacitance  (bottom  right)  the  waveforms  become  mirror-symmetric. 
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Figure  4.10:  Selected  simulations  from  previous  figure,  comparing  tip  voltage  output  be¬ 
fore  and  after  processing  through  an  output  non-linearity  (sigmoid  shown  in  Rgure  4.12). 
Reversal  of  OS  is  seen  for  the  simulation  on  the  second  row,  in  which  the  sigmoid  gen¬ 
erates  a  larger  total  output  for  the  “fatter”  leftward  motion  response.  Under  the  same 
conditions,  but  with  negligible  capacitance  (bottom  row),  the  sigmoid  merely  preserves 
the  mirror-symmetry  of  the  voltage  output. 
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Figure  4.11:  DI  for  the  linear  average  of  distal  tip  waveforms  obtained  from  simulations 
of  cell  in  Figure  4.8,  over  a  range  of  inhibitory  synaptic  conductance  densities  and 
stimulus  velocities,  for  Cm  =  1.0/iF/cxn^  (middle  left)  and  10~‘fiF/cm’  (lower  right). 
Axes  are  shown  in  the  upper  right.  For  this  linear  functional  DS  is  eliminated  as  inhibition 
is  lowered,  independent  of  Cm.  Furthermore,  cable  capacitance  does  not  seem  to  have 
a  significant  effect  on  DS  when  inhibition  is  present. 


Figure  4.12:  DI  for  the  same  simulations  as  Figure  4.11,  but  taken  of  the  average  of 
distal  tip  waveforms  after  being  passed  through  the  sigmoid  shown  in  upper  right.  With 
normal  Cm  (middle  left)  the  directionality  reverses  {DI  <  0)  as  inhibition  is  lowered. 
When  Cm  and  inhibition  are  both  small  (lower  right)  this  reversal  is  dramatically  reduced. 
Axes  are  as  Figure  4.11.  Although  the  strength  of  the  DS  reversal  is  rather  small  for  the 
particular  parameters  in  these  simulations,  subsequent  non-linearities  could  in  principal 
amplify  the  predicted  negative  DI. 


Figure  4.13:  Geometry  of  apparent  motion  simulations  using  the  starburst  amacrine  celt 
described  earlier,  with  the  addition  of  an  inactivating  A"*'  channel  at  each  of  the  distal 
portions  of  the  tree.  Inset  shows  the  time  course  of  the  preferred  direction  apparent 
motion  sequence.  As  before,  the  tip  output  of  interest  is  the  highlighted  node  at  the 
right  edge. 


Figure  4.14:  Voltage  at  the  righthandmost  distal  tip  in  Figure  4.13,  in  response  to 
slit  A  and  B  alone  and  in  the  preferred  and  null  direction  apparent  motion  sequences. 
Preferred  direction  facilitation,  in  which  the  amplitude  of  the  response  to  B  is  enhanced 
by  the  prior  stimulus  A,  is  due  to  the  voltage-dependent  inactivation  of  a  distal 
channel  by  A.  Null  direction  inhibition,  in  which  the  amplitude  of  the  response  to  A  is 
reduced  by  the  prior  stimulus  B,  is  due  to  the  on-the-path  shunting  described  in  the 
text. 
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Chapter  5 

Discussion 


We  shall  now  recapitulate  and  discuss  the  experimental  and  theoretical  re¬ 
sults,  and  discuss  some  related  issues,  including  that  of  developmental  con¬ 
siderations  and  the  possible  relevance  of  the  model  and  data  to  other  brain 
areas  and  computations. 


5.1  Experimental  Evidence  of  Pre- Synaptic  Com¬ 
putation  of  Directional  Selectivity  and  Pos¬ 
sible  Complications 

Again,  we  interpret  the  recordings  shown  in  Chapter  3  as  evidence  for  ex¬ 
citatory  synaptic  input  onto  DS  ganglion  cells  in  turtle  retina  as  being  al¬ 
ready  DS.  This  finding  pushes  the  computational  elements  for  DS  back  onto 
amacrine  cells  or  further.  Thus,  in  Chapter  4  we  have  postulated  a  specific 
model  for  DS  in  which  evenly  distributed  excitatory  and  (slower)  inhibitory 
inputs  onto  individual  amacrine  cell  dendrites  generate  DS  excitatory  out¬ 
puts  onto  DS  ganglion  cells,  the  outputs  being  at  or  near  the  tips  of  the 
amacrine  cell  dendrites.  The  crucial  asymmetry  in  this  model  is  found  in 
the  receptive  fields  (the  spatial  extent  of  the  input)  of  the  distal  tip  output 
synapses:  the  input  onto  the  amacrine  cell,  although  evenly-distributed  and 
with  symmetric  receptive  fields,  is  displaced  with  respect  to  any  tip  out¬ 
put.  This  asjrmmetry  must  be  maintained  in  the  aiferents  onto  any  given 
DS  ganglion  cell,  that  is  the  ganglion  cell  must  contact  amacrine  cell  tips 
of  a  specific  orientation  (see  Section  5.5).  Also,  the  most  important  non¬ 
linearity  in  the  system  arises  &om  the  interaction  between  the  excitatory 
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and  inhibitory  synaptic  conductance  input  onto  the  amacrine  cell  dendrites. 

The  observation  that  under  our  recording  conditions  IPSPs  are  reduced 
soon  after  whole-cell  access  supports  the  above  interpretations,  but  we  may 
ask  whether  this  blockage  is  complete  for  distal  inhibitory  inputs  and,  if 
not,  whether  such  inputs  may  preserve  DS  via  a  correlative  null  direction 
vetoing  mechanism.  If  a  synaptic  event  (such  as  distal  inhibition)  is  able 
to  influence  somatic  events,  however,  then  in  principle  that  event  must  be 
observable  by  somatic  recordings  (though  not  necessarily  accxirately  mea¬ 
sured).  Since  neither  inhibitory  voltages  (under  current  clamp)  nor  currents 
(imder  voltage  clamp)  were  observed  in  the  motion  responses,  then  we  sub¬ 
mit  that  the  blocking  of  inhibitory  input  in  our  experiments  is  sufficient  to 
support  the  DS  excitatory  input  interpretation. 

A  similar  argument  may  be  used  if  we  consider  possible  contributions 
from  other  DS  ganglion  cells  which  do  not  have  their  inhibition  blocked  but 
eire  able  to  transmit  DS  potentials  via  gap  junctions  to  the  recorded  cell.  In 
this  case,  as  well,  the  electrical  continuity  of  the  gap  junctions  would  mean 
that  any  inhibitory  event  that  could  eventually  influence  the  soma  voltage 
of  the  recorded  cell  must,  in  turn,  be  observable  firom  that  cell. 

5.2  Anatomical  Predictions  of  the  Cable  DS  Model 

The  essence  of  the  DS  model  presented  'n  Chapter  4  is  that  amacrine  ceU 
output  synapses  that  are  asymmetric  with  respect  to  the  dendritic  tree  of  the 
amacrine  cell  will  have  directional  responses  which,  under  reasonable  condi¬ 
tions,  will  be  DS.  We  have  focused  on  the  distal  tip  sjmapses  in  particular 
since,  assuming  that  the  input  to  the  tree  is  homogeneous,  these  locations 
have  the  most  asymmetric  input.  Since  a  given  DS  ganglion  cell  wotild  col¬ 
lect  DS  inputs  of  a  given  orientation,  then  the  model  predicts  an  overall 
spatial  asymmetry  in  the  amacrine  cell  connectivity  onto  the  ganglion  cell. 
For  the  model  version  discussed  at  length,  the  input  cells  would  be  shifted 
towards  the  null  side  of  the  ganglion  cell  receptive  field  (see  Figure  5.1).  It 
may  also  be  that  the  amacrine  ceU  DS  outputs  pass  through  bipolar  cells 
before  going  to  the  ganglion  cell:  for  the  model  version  discussed  here  the 
amacrine  cell  asymmetry  is  still  towards  the  null  side  for  a  sign-preserving 
bipolar  connection,  and  towards  the  preferred  side  for  a  sign-inverting  bipo¬ 
lar  connection.  DS  facilitatory  data  (Section  4.1.5),  however,  argue  agiunst 
the  latter  possibility. 
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Preferred  Direction 


Figure  5.1:  The  model  of  DS  amacrine  cell  tip  excitatory  outputs  predicts  that  the  den¬ 
dritic  fields  of  amacrine  cells  presynaptic  to  a  given  DS  ganglion  cell  would  be  displaced 
towards  the  null  direction,  with  respect  to  the  ganglion  cell's  dendritic  field  (see  also 
(661). 
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5.2.1  Correlating  Anatomy  with  Function 

One  experiment  to  test  these  predictions  wotild  be  to  inject  a  physiologically 
characteri2ed  DS  ganglion  cell  with  a  retrograde  tracer,  such  eis  Herpes  Sim¬ 
plex  virus  [52],  and  then  compare  the  distribution  and  orientation  of  stained 
dendritic  afferents  with  the  DS  receptive  field  properties.  The  whole-cell 
technique  described  here  may  be  well  suited  for  this  type  of  experiment  be¬ 
cause  a)  the  wide  bore  of  the  patch  electrode  facilitates  delivery  of  electrode 
contents  to  the  cell  and  h)  the  stability  of  the  method  allows  for  detailed 
mapping  of  the  receptive  field. 


5.3  Relative  Frequency  of  Directional  Selective 
Responses 

At  this  stage  of  otir  work  we  have  avoided  characterizing  the  percentage 
of  DS  cells  using  this  preparation,  primarily  because  a)  the  protocols  and 
recording  techniques  were  evolving,  b)  the  DS  criteria  described  in  Chapter  3 
are  probably  weaker  than  that  usually  used  (explicitly  or  implicitly)  in  the 
literature,  and  c)  for  many  cells  a  DS  response  from  a  WF  protocol  was 
not  stable  over  time.  With  respect  to  the  second  point,  while  our  criteria 
does  pick  out  statistically  significant  asymmetries,  it  does  not  consider  the 
strength  of  the  as3rmmetry  (e.g.  DI).  With  these  pointr  m  r’^nd,  for  the 
cohort  oi  on  ceils  described  in  Appendix  B,  si  50%  showed  DS  responses  at 
some  point  during  the  recording.  The  cells  shown  in  Chapter  3  were  chosen 
because  they  had  particularly  stable  DS  responses.  As  we  have  seen,  all  of 
these  cells  show  evidence  for  DS  synaptic  input.  On  the  other  hand,  we 
have  found  no  cells  in  which  an  initial  DS  response  disappeared  over  time 
(i.e.  subsequent  to  removal  of  inhibition  at  the  cell)  without  a  concurrent 
degradation  of  the  cell  response  in  general.  In  other  words,  we  found  no 
evidence  that  ganglion  cell  DS  in  turtle  depends  on  the  presence  of  inhibition 
onto  that  ceU. 
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5.4  Comparison  to  Results  of  Watanabe  and  Mu¬ 
rakami:  Directionally  Selective  Responses  in 
ON,  OFF,  and  ON/OFF  Cells 

Our  experimental  results  are  similar  to  that  of  Watanabe  and  Murakami 
([67]).  However,  in  addition  we  have  observed  evidence  for  DS  input  into 
DS  ON/OFF  ganglion  cells,  as  well  as  into  ON  and  OFF  (not  shown)  DS 
ganglion  cells.  These  results  do  not  necessarily  contradict  Watanabe  and 
Murakamis’  observations,  since  in  their  paper  they  simply  did  not  observe 
IPSPs  at  all  in  either  ON  nor  OFF  ceUs.  It  is  quite  possible  that  if  they 
were  either  able  to  selectively  block  inhibitory  input  to  the  ganglion  cell, 
or  if  they  voltage  clamped  their  cells,  as  we  have  done  here,  then  they  may 
have  concluded  that  the  observed  IPSPs  for  motion  responses  in  ON/OFF 
cells  were  not  necessary  for  DS. 


5.5  Development  of  Direction  Selectivity:  The 
Problem  of  Coordination  of  Asymmetries 

A  problem  for  any  model  of  DS  is  how  to  break  symmetry  on  a  scale  sig¬ 
nificantly  larger  than  that  predicted  by  random  distributions  of  local  asym¬ 
metries.  In  the  r'tinn,  the  DS  response  is  correlated  over  multiple  subunits, 
and  several  directions  are  represented,  hor  the  model  presented  here,  this 
translates  into  determining  how  a  ganglion  cell  could  selectively  connect 
to  amacrine  dendrite  tips  of  similar  orientation.  One  solution  that  seems 
natural  for  this  model  is  to  postulate  a  Hebbian-type  correlational  process 
([35])  similar  to  long-term  potentiation  in  the  hippocampus  (see  review  by 
Brown,  Kairiss,  and  Keenan  [15]).  A  Hebbian-type  process  strengthens  ac¬ 
tive  synapses  when  the  pre  and  postsynaptic  cells  fire  simultaneously  and 
weakens  synapses  when  activity  is  not  correlated  (see  Figure  5.2). 

The  initial  syxnmetry  in  the  retina  may  then  be  broken  as  follows:  As¬ 
sume  that  at  an  early  stage  of  synaptogenesis  initial  DS  amacrine  connec¬ 
tions  to  a  proto-DS  ganglion  cell  have  a  range  of  orientations,  but  overall 
there  is  a  slight  orientation  bias.  If  this  bias  is  strong  enough,  then  a  Heb- 
bian  mechanism  could  reinforce  it.  The  connectivity  of  subsequent  synapses 
would  then  become  stronger  or  weaker  depending  on  their  dendrite’s  align¬ 
ment  relative  to  the  initial  (weak)  orientation. 

A  Hebbian  mechanism,  at  least  in  this  simple  form,  is  suited  for  a  pre- 
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gangUonic  DS  model  because  the  mechanism  lequixes  signals  -which  can  be 
correlated;  the  postsynaptic  cell  tends  to  collect  correlatable  inputs.  Gan¬ 
glionic  computation  of  DS  loses  this  advantage,  since  there  is  no  motion 
asymmetry  intrinsic  in  the  ganglion  cell  aiferents. 

This  scheme  requires  the  proper  IPL  input,  that  is  motion  stimuli,  to 
work.  Thus  it  will  be  important  to  determine  when  DS  properties  emerge 
relative  to  visual  experience  in  any  given  species.  It  is  possible,  though, 
that  prenatal  DS  properties  may  be  formed  in  response  to  motion~$imidating 
stimuli  onto  the  IPL,  for  example  spontaneous  retinal  waves  [51]  (see  Section  5.7). 


5.6  Retinal  Directional  Selectivity:  Exemplar  of 
a  Canonical  Computational  Mechanism? 

Thus  far  we  have  discussed  the  directional  properties  that  may  arise  intrin¬ 
sically  whenever  dendrites  have  distributed  excitatory  and  inhibitory  input, 
at  least  in  terms  of  the  waveforms  at  the  dendrite  end.  The  key  to  this  prop¬ 
erty  is  the  asymmetric  distribution  of  the  inputs  with  respect  to  the  output, 
and  we  have  demonstrated  how  this  sort  of  asymmetry  may  be  especially 
effective  when  the  output  synapse  is  on  the  dendrite  tip.  Since  in  retina  the 
salient  stimuli  features  are  retinotopic,  it  is  very  important  to  consider  the 
exact  geometry  of  the  amacrine  cells’  dendritic  trees. 

5.6.1  For  Generic  Central  Neurons:  Trees  of  Dendrite  Filters 

However,  in  the  more  general  case,  we  may  consider  the  cascade  of  inputs 
along  a  dendritic  branch  of  a  generic  central  neuron.  The  performance  of  the 
retinal  DS  model  suggests  that  a  similar  directionality  will  exist  under  some 
conditions  for  non-retinal  neurons,  with  respect  to  the  aiferents  along  each 
branch.  We  suggest  that  with  plausible  biophysical  parameters  the  soma  of 
a  more  general  cell,  which  is  asymmetric  with  respect  to  the  dendritic  input, 
will  see  a  ‘‘directional”  response  from  each  of  its  branch  inputs.  Specifically 
(at  least  with  the  version  of  the  model  disctissed  here),  for  a  given  set  of 
inputs  over  time  along  a  branch,  the  soma  EPSP  will  be  greatest  when  the 
temporal  order  of  the  inputs  is  from  distal  to  proximal,  i.e.  stepping  in  time 
toward  the  soma.  Conversely,  the  same  set  of  inputs,  only  reversed  in  time, 
may  be  much  less  effective  in  depolarizing  the  soma. 

With  the  conventional  neuron’s  dendrite-to-soma  signal  flow,  the  polarity 
of  the  filter  is  reversed  from  that  of  the  retinal  version:  the  relevant  output 
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Random  Connections  with 
Slight  Overall  Bias 
Pointing  Rightwards 


»  I  \ 

*  •  > 


Connection  Bias  Reinforced 
by  Hebbian  Mechanism 


Figrire  5.2:  According  to  the  model  presented  here,  Hebbian  mechanisms  may  be 
used  during  development  to  reinforce  an  initial  slight  directional  bias  between  random 
amacrine  ceil  DS  tip  connections  and  a  proto-DS  ganglion  cell.  This  DS  tip  output 
model  seems  well  suited  for  such  a  scheme  since  there  is  a  directly  correlatable  compo¬ 
nent  between  the  pre-synaptic  signals.  In  the  mature  circuit  on  the  right,  the  dashed 
amacrine  ceils  have  all  of  their  connections  weakened. 


of  a  given  branch  is  its  proxinoal  (somatic)  end.  However,  as  with  retinal 
directioned  selectivity,  the  net  result  is  that  each  dendritic  branch  functions 
as  more  than  just  a  time-independent  “point”  integrator  of  inputs.  Rather, 
the  bramches  fimction  as  nonlinear  “spatio”-temporal  filters.  We  may  also 
consider  implications  of  trees  of  such  filters,  with  respect  to  the  somatic 
response  to  the  cascading  of  branch  filter  characteristics.  Finally,  we  note 
that  this  model  allows  for  the  efficient  implementation  of  various  filters  which 
operate  on  the  same  input  “bus”,  with  the  preferred  sequence  of  each  filter 
determined  solely  by  the  postsynaptic  organization  (see  Figure  5.3). 


Figure  5.3:  An  afferent  fiber  tract  could  form  synapses  onto  dendrites  of  different  cells 
with  various  spatial  ordering,  relative  to  the  soma  of  each  cell,  thus  forming  various 
non-linear  spatio-temporal  filters.  In  this  cartoon  the  three  target  neurons  would  have 
different  preferred  temporal  sequences  of  the  signals  A,  B  and  C  off  the  input  bus. 


In  the  more  general  case  the  coding  of  the  inputs  onto  a  dendrite  depends, 
of  course,  on  the  cell.  One  possible  candidate  would  be  phonemes  from  an 
auditory  input  stream  in  which  the  temporal  order  of  the  phonems  is  unique 
for  a  given  word. 


98 


5.7  So  What  Are  Ganglion  Cells  For,  Anyway? 

The  resiilts  of  the  work  presented  here  suggest  that  the  DS  ganglion  cell  in 
turtle  doesn’t  do  very  much  in  terms  of  the  DS  computation  itself;  rather, 
the  ganglion  cell’s  job  seems  to  be  to  collect  many  instances  of  a  specific 
computation  -  done  by  amacrine  cells  -  in  order  to  form  a  coherent,  consistent 
receptive  field.  It  may  be  that  this  sort  of  final  “editing”  of  retinal  processing 
is  typical  for  the  ganglion  cell  layer. 

As  we  mentioned  above,  the  conditions  for  DS  dendritic  waveforms  are 
weaik:  therefore  we  suggest  that  not  only  may  this  mechanism  be  relevant 
outside  the  retina,  but  also  that  some  degree  of  DS  may  be  the  rule  for 
IPL  output  rather  than  the  exception.  As  we  mentioned  in  Section  4.1.1, 
if  we  assume  that  amacrine  cells  in  general  have  outputs  throughout  their 
dendritic  trees,  then  most  of  the  outputs  will  be  asymmetric  with  respect  to 
the  inputs  that  can  effect  them.  The  degree  of  input  asymmetry  for  a  ^ven 
output  may  in  turn  be  reflected  in  the  degree  of  DS  for  that  output.  If  this  is 
the  case,  then  we  may  surmise  the  following:  a  given  species  will  demonstrate 
retinal  DS  if  either  a)  the  ganglion  cells  have  the  necessary  machinery  (e.g. 
Hebbian)  to  selectively  collect  DS  IPL  signals  and/or  b)  motion  or  motion- 
simulating  input  onto  the  IPL  is  present  during  a  critical  plastic  stage  of 
development.  This  leads  to  the  hypothesis  that  the  existence  of  retinal 
DS  for  a  species  depends  more  on  detaUs  of  developmental  mechanisms  and 
timing  fc-  that  species,  rather  than  on  any  “special”  biophysical  m«;vua.....uj 
involved  in  the  mature  DS  circuitry. 


5.8  Implications  for  DS  in  the  Visual  Systems  of 
Higher  Species 

As  we  mentioned  in  Chapter  2,  DS  responses  in  higher  species  tend  to  first 
emerge  in  visual  cortex.  Does  the  model  presented  here  have  any  relevance 
for  cortical  DS?  Certainly  visual  cortex  has  the  retinotopic  organization  that 
makes  this  model  arise  in  a  simple  way  within  retinal  tissue.  Similar  bio¬ 
physical  mechanisms  probably  exist  in  cortex  as  well,  as  discussed  above. 
Whether  or  not  the  S3maptic  organization  necessary  to  recapitulate  the  reti¬ 
nal  DS  model  exists  in  cortex  remiuns  to  seen,  but  hopefully  this  model 
will  generate  testable  predictions  for  DS-generating  connectivity  in  visual 
cortex. 
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Appendix  A 

Isolation  and  Recording  of 
Turtle  Retina  Allowing  for 
Arbitrary  Electrode 
Approach  and 

Trans-Illumination  of  Retina 


Isolated  retina  pre*  rations  have  proved  useful  in  many  experimental  pro¬ 
tocols  since  developed  by  Ames  and  Nesbett  for  rabbit  retina  ([!])■  These 
advantages  include: 

•  Better  clearance  properties  (useful  in  pharmacological  studies). 

•  Avoidance  of  electrode  shadowing  of  stimulus. 

•  Infra-red  trans-illumination  of  the  tissue  which  facilitates  electrode 
placement  and  control. 

•  Visual  observation  during  recording  which  allows  for  impalement  of 
selected  (e.g.  stained)  cells. 

In  1989  Lipets  et.  al.  [45]  described  a  new  technique  for  the  isolation  of 
a  functional  turtle  retina.  While  their  technique  apparently  yields  normal 
responses  from  cones  and  horizontal  cells,  it  does  not  allow  either  direct 
electrode  approach  from  the  ganglion  cell  side,  nor  transmitted  light  mi¬ 
croscopy.  In  this  appendix  we  describe  a  modification  of  their  technique, 
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which  2Jlows  both  physical  and  visual  access  to  either  side  of  the  retina 
([11]).  In  particular,  in  the  majority  of  the  dissections  it  is  possible  to  lo¬ 
cate  the  visual  streak  and  use  this  as  a  location  reference  throughout  the 
experiment.  Because  we  initially  used  extracellular  recordings  to  assess  this 
preparation  (with  metal  electrodes),  and  because  intracellular  protocols  will 
be  discussed  in  Appendix  B,  in  this  appendix  we  will  confine  recording  re¬ 
sults  to  conventional  extracellular  protocols. 


A.l  Dissection  Method 

Prior  to  the  dissection  a  small  piece  of  filter  paper  is  prepared  and  mounted 
in  a  holder  (Figures  A.l  and  A.2).  The  arm  of  the  holder  is  mounted  on  the 
head  carriage  of  a  dissecting  microscope  (with  the  optics  removed),  which 
serves  as  a  useful  platform  for  the  dissection.  After  mounting,  the  paper 
holder  is  adjusted  so  that  the  filter  paper  is  coplanar  with  the  microscope 
stage. 

With  respect  to  the  filter  paper,  we  have  had  some  success  with  both 
coffee  maker  filter  paper  (Mellita)  and  nitro-ceUulose  paper  (Millipore,  0.45 
fiva.  pore  size).  Generally,  the  nitro- cellulose  paper  seems  to  have  slightly 
more  affinity  for  the  retina,  and  it  is  a  simple  matter  to  melt  small  windows  in 
this  material  with  the  heated  end  of  glass  tubing  of  the  appropriate  diameter. 

Perfusate  is  taken  from  [5]:  ST. 5  mM  NaCl,  2.6  mM  KCl,  2  mM  MgCljt 
4  mM  CaClj,  31.5  mM  NaHCOs,  10  mM  glucose,  10  mM  HEPES,  bubbled 
with  95%03/5%C03.  In  practice,  10 X  stock  solutions,  without  CaCl3  or 
glucose,  were  prepared  ahead  of  time  and  stored  in  the  refrigerator. 

Turtles  are  dark-sulapted  in  a  light-tight  box  for  at  least  8  hours  prior  to 
an  experiment.  After  adaptation  the  entire  dissection  procedure,  including 
decapitation,  is  done  using  a  head-mounted  monocular  infra-red  (IR)  viewer 
(Model  6100M,  Electrophysics  Corporation,  Nutley,  New  Jersey).  After  de¬ 
capitation  and  pithing  the  head  is  immediately  mounted  in  a  dissection 
holder,  covered  with  foil  to  block  out  light,  and  put  in  the  refrigerator  (S^C). 

After  a  variable  amount  of  time  (ranging  from  immediately  to  one  day) 
the  head  and  holder  are  unwrapped  and  one  eye  is  enucleated.  The  head 
is  then  wri^ped  back  in  foil  and  returned  to  the  refrigerator  for  possible 
dissection  of  the  remaining  eye  later.  The  enucleated  eye  is  placed  on  blotter 
paper  and  hemisected,  and  the  posterior  half  laid  face  down  on  the  paper 
to  drain  for  about  20  seconds.  The  drained  eyecup  is  then  rinsed  in  the 
Ringer’s  solution  for  about  one  minute.  The  eyecup  is  then  removed  from 
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the  Ringer’s  Jind  the  disc  and  streak  are  located.  The  cup  is  bisected  on 
a  non-porous  surface  (to  avoid  premature  removal  of  the  retina)  along  the 
perpendicular  to  the  streak,  narrowly  skirting  the  disc.  One  half  of  the 
bisected  cup  is  placed  in  a  small  dish  of  the  Ringer’s  solution  and  put  into 
the  refrigerator  for  later  use  (up  to  8  hours). 

If  the  remaining  half  has  the  disk  then  another  cut  is  made,  parallel  to 
the  first,  to  remove  the  disc  from  the  section.  Further  trimming  may  be 
done  in  order  to  make  the  section  about  4mm  wide.  If  necessary,  the  section 
is  Ladd  on  its  side  on  the  non-porous  surface  in  order  to  dradn  more  vitreous. 
The  section  is  then  plau:ed  on  blotter  paper,  amd  a  scalpel  is  used  to  make 
two  lateral  cuts  (perpendictilar  to  the  long  auds)  through  the  retina  only, 
dividing  the  section  in  thirds.  The  section  is  then  centered,  retina  side  up, 
on  a  small  (ss  3mm  x  5mm  x  1mm  high,  made  of  tape)  support  which  is 
mounted  on  the  XY  stage  of  a  dissection  microscope.  The  laterad  ‘‘ears” 
of  the  section  are  held  down  by  two  spring  wire  clips  such  that  the  middle 
third  of  the  retina  is  elevated  slightly  by  the  support. 

The  filter  paper  retina  catrrier  hole  is  then  adigned  in  the  XY  plame,  with 
the  paper  several  mm  above  the  tissue,  so  that  the  hole  is  centered  on  the 
middle  third  of  the  section.  Monitoring  from  the  side,  the  paper  is  then 
lowered  using  the  dissection  microscope’s  focus  adjustment  tmtil  the  paper 
bairely  rests  on  the  retina.  The  paper  is  then  adlowed  to  sit  on  the  retina  for 
about  20  to  40  seconds,  and  then  the  paper  '  'der  is  lifted  quickly  by  hand, 
taddng  the  section  of  retina  with  it.  The  retina  section  is  now  attached  by 
its  gamglion  cell  surface  to  the  underside  of  the  filter  paper.  At  this  point 
the  dissection  is  inspected  and,  if  it  is  not  satisfactory,  adjusted  with  fine 
forceps  in  room  light.  In  our  experience,  light  exposure  does  no  harm  once 
the  retina  is  isolated.  The  retina  and  filter  paper  is  then  placed  under  a 
nylon  ring  in  the  perfusion  bath  (flow  «lml/min).  Recordings  are  begun 
after  at  least  a  one  and  half  hour  wait,  but  consistent  light  responses  have 
been  recorded  within  minutes  of  the  dissection. 

A. 1.1  Importance  of  Low  Light  During  Dissection 

The  key  step  to  this  technique  is  that  of  maintaining  very  low  light  levels 
during  the  dissection  to  prevent  the  interdigitation  of  the  pigment  epithelium 
with  the  photoreceptors  prior  to  and  during  the  actual  isolation  [17].  Other 
factors  that  may  contribute  to  retinal  adhesion  to  the  pigment  epithelium, 
including  temperature  or  pH  [70],  were  not  evaluated.  We  have  fotmd  that 
refrigeration  schedules  are  not  critical,  and  the  material  is  relatively  resistant 
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to  2uioxia,  since  we  have  obtained  apparently  reasonable  light  responses  from 
both  freshly  killed  emimals  and  from  heads  that  have  been  stored  up  to  two 
full  days  under  refrigeration. 


A.2  Recording  System  for  Superfusion  and  Transmitted 
Light  IR  Microscopy 

The  perfusion  system  and  stimulus/recording  layout  is  set  up  to  provide 
mixed  stimulus  and  IR  illumination  from  below,  and  recording  and  IR  video 
monitoring  from  above,  as  shown  in  Figure  A.3. 


A.3  Stimulus  Generation,  Data  Acquisition  and 
Experiment  Control 

Light  stimuli  were  generated  on  a  Teztronix  608  monitor  with  a  P4  phos¬ 
phor,  controlled  by  a  Picasso  CRT  Image  Synthesizer  (Innisfree,  Cambridge, 
Massachusetts).  All  experiments  described  in  this  thesis  used  mavirrmtin  con¬ 
trast  stimuli  (full  brightness  firom  full  off).  The  monitor  image  was  collected 
by  a  4  inch  diameter  convex  lens  and  then  mixed  into  the  microscope  IR 
illuminator  path  by  a  half-silvered  mirror  prior  to  the  microscope  condenser. 
The  final  image  on  the  retina  was  demagniiied  by  a  fac  er  of  40,  and  the 
intensity  at  retina  was  »  25/iW/cm^  (estimate). 

Control  of  the  Picasso  and  data  acqtiisition  was  handled  with  an  IBM 
Model  AT  Personal  Computer,  running  under  the  Asyst  system  (Macmillan 
Software  Company,  New  York,  New  York).  Electrophysiological  amplifier, 
window  discriminator  (World  Precision  Instruments),  and  Picasso  phase  out¬ 
puts  were  both  recorded  on  a  PCM  VCR  recorder  (Vetter),  and  digitized 
for  on-line  data  analysis  and  storage.  Off  line  data  analysis  was  done  with 
various  standard  computing  resources. 


A.4  Extracellular  Recordings 

In  order  to  assess  the  viabihty  of  this  preparation,  extracellular  recordings 
were  done,  presumably  picking  up  spiking  activity  from  ganglion  cells  (Da- 
gan  extracellular  amplifier,  Dagan  Corporation,  Minneapolis,  Minnesota). 
Cells  were  located  by  advancing  extracellular  electrodes  (2-5Mn,  parylene 
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coated  tungsten  electrodes,  A-M  Systems,  Everett,  Washington)  under  vi¬ 
sual  control  onto  the  inner  surface  of  the  retina.  Extracellular  spikes  (typ¬ 
ically  50  -  200  /iV)  were  obtained  both  in  response  to  flashing  spot  stimuli 
(typically  400  /im  diameter,  centered  on  the  electrode,  Ihx)  and  sponta¬ 
neously.  Depending  on  the  experiment,  spikes  were  found  both  when  the 
electrode  tip  contacted  the  surface  of  the  retina,  and  sometimes  only  when 
the  electrode  was  advanced  (at  a  45°  angle  to  the  surface)  %  50  —  200/im 
beneath  the  surface.  In  the  latter  case,  severe  indentation  of  the  surface  was 
seen. 


A.5  Properties  of  Ganglion  Cells 

A. 5.1  ON,  OFF,  and  ON/OFF  Properties 

Basic  ON,  OFF,  and  ON/OFF  ganglion  cell  physiological  types  have  been 
recorded  in  this  preparation.  Results  from  extracellular  experiments  on  a 
group  of  45  cells  are  shown  in  Table  A.l.  In  these  cells,  the  receptive  field 
was  mapped  with  a  Ihz  200  by  200  tiro,  spot,  stepping  over  a  900  by  900  fim. 
grid,  and  were  classified  according  to  the  dominant  response.  Some  cells 
showed  a  prominent  ON  or  OFF  response  in  the  majority  of  their  spot 
receptive  field,  but  would  show  a  clear  reversal  of  type  in  restricted  areas. 
These  cells  are  listed  in  the  table  as  ** Variable”. 


Type 

Number 

Percentage 

ON 

8 

OFF 

21 

ON/OFF 

13 

VARIABLE 

8 

Table  A.l:  Basic  Ganglion  Cell  Types  from  Extracellular  Recordings  in  the 
Isolated  Intact  Turtle  Retina  (N  =  45) 


A.5.2  Receptive  Field  Size  and  Location,  Center  Surround 
Organization 

The  distribution  of  receptive  field  si2es  of  the  cohort  of  cells  in  Table  A.l 
is  illustrated  in  Figure  A.4.  Field  properties  were  mapped  out  as  described 
above,  and  a  Student  T  test  was  applied  to  the  response  at  each  grid  point, 
selecting  responses  that  were  significantly  above  background  at  the  5  percent 
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confidence  level.  The  total  area  of  the  grid  covered  by  the  selected  points 
was  then  calculated,  and  the  diameter  of  a  circular  field  with  the  same  area 
was  derived. 

Visuedization  of  the  retina  during  recording  allows  for  immediate  de¬ 
termination  of  the  relation  between  soma  location  and  the  receptive  field, 
relative  to  the  visual  streak.  In  our  experiments  we  have  consistently  ob¬ 
served  a  bias  in  this  relation  which  tends  to  put  the  receptive  field  center 
midway  between  the  electrode  tip  and  the  streak. 

Some  cells  showed  a  clear  excitatory-center,  inhibitory-surrouiid  response 
to  spots  of  increasing  diameter.  These  characteristics  aie  discussed  in  detail 
in  Appendix  B. 


A. 6  Preservation  of  Photoreceptors 

A  major  concern  with  the  isolation  procedure  is  the  integrity  of  the  frag¬ 
ile  photoreceptor  outer  segments.  Scaxming  electron  micrographs  of  the 
photoreceptor  surface  is  one  way  to  evaluate  this  morphology  [45].  In  the 
micrographs  shown  in  Figiixe  A.5  ‘•he  outer  segments  ate  clearly  visible,  and 
large  areas  of  photoreceptors  seem  to  have  weathered  the  dissection  rather 
well.  We  also  saw  extensive  areas  of  damage,  but  we  assume  that  the  “good” 
areas  are  large  enough  for  the  electrophysiology. 


A.7  Infra-Red  Viewing  of  Preparation  During  Record¬ 
ing 

An  important  advantage  of  this  preparation  is  the  ability  to  see  the  tissue 
during  recording,  especially  for  the  whole-cell  patch  protocol  described  in 
the  following  appendix.  An  example  IR  video  hardcopy  of  a  retina  during  a 
whole-cell  protocol  is  shown  in  Figure  A.6. 


A. 8  Consistency  of  Method 

This  method  has  yielded  preparations  in  which  extracellular  ganglion 
responses  are  quite  readily  obtained.  In  addition,  there  appears  to  be  no 
obvious  large-scale  artifact  in  the  mapped  receptive  fields  that  might  be  due 
to  the  dissection. 
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A. 9  Stability  and  Longevity  of  Method 

We  have  not  encountered  significant  bleaching  effects  with  the  isolated  prepa¬ 
ration,  despite  the  absence  of  pigment  epithelium,  possibly  due  to  the  low 
level  of  illumination.  Vigorous  responses  to  spots  have  been  obt^dned  from 
retin2is  after  up  to  14  hours  of  a  protocol  consisting  of  200/im  diameter  flash¬ 
ing  spots  or  moving  gratings,  presented  at  various  locations  over  a  2  x  2mm 
retinal  area  for  %  80%  of  the  time,  and  large  (500  —  lOOO^m  diameter)  stim¬ 
uli  for  the  remaining  time.  In  addition,  extracellular  recordings  from  single 
cells  have  been  maintained  for  up  to  8  hours.  This  stability  is  especially  use¬ 
ful  when  investigating  the  Ane  structure  of  the  terrapin  retina’s  complicated 
receptive  fields. 

A. 10  Summary 

We  have  demonstrated  that  fully  isolated  turtle  retina  is  a  viable  prepara¬ 
tion.  The  dissection  technique  that  we  describe  is  designed  to  be  as  robust 
as  possible,  and  we  believe  that  this  preparation  is  amenable  to  a  variety  of 
extracellular,  intracellular,  pharmacological,  and  histological  techniques. 
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Figure  A.l;  Holder  for  filter  paper  retina  carrier  allows  alignment  of  paper  with 
dissecting  microscope  sia^vj  prior  to  dissection  under  IR.  Modified  alligator  clip 
holds  paper.  Holder  arm  is  attached  to  head  of  microscope  for  precision  control 
of  Z-axis  movement. 
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Figure  A.2:  Procedure  for  the  dissection  of  intact  isolated  turtle  retina.  Inverted 
retina  section  is  mounted  on  dissecting  microscope  stage  for  precision  control 
of  XY-axis  movement.  Microscope  head  optics  are  removed  and  visualization 
from  above  (for  aligning  paper  with  tissue)  and  from  the  side  (for  lowering  of 
paper  onto  tissue)  is  only  through  IR  goggle.  Detached  retina  and  paper  are 
immediately  mounted  under  nylon  ring  in  perfusion  bath  (See  Figure  A.3). 
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Figure  A.3:  Layout  of  tissue,  bath,  stimulus  optics,  infra-red  illumination,  and 
electrode  during  recordings  from  intact  isolated  retina. 
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Figure  A.4;  Distribution  of  diameters  of  circles  with  equivalent  areas  of  recep¬ 
tive  fields  for  cells  described  in  Table  A.l  {N  =  45). 


Figure  A.5:  Scanning  electron  micrographs  of  photoreceptor  side  of  isolated  intact 
turtle  retina.  Top:  Large  areas  of  photoreceptors  do  not  show  obvious  damage  after  the 
dissection  described  in  this  appendix  (magnification  «  350a,  zoom  box  %  80  x  50/im). 
Bottom:  Blow-up  of  the  zoom  box  in  top  figure  (magnification  st  1500a),  showing  a 
tangle  of  intact  photoreceptor  outer  segments. 
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Figure  A.6:  Retina  as  seen  during  recording  (Cell  A67).  This  scanned  image  was  taken 
from  a  video  hardcopy  of  the  IR  video  monitor  (final  magnification  80z).  The  whole-cell 
patch  electrode  enters  from  the  upper  left  corner,  and  the  visual  streak  courses  vertically 
on  the  left  hand  side,  coming  to  within  250/im  of  the  electrode  tip.  The  receptive  field 
of  this  cell  was  %  300/im  in  diameter,  and  was  centered  lOO/im  to  the  left  and  200/im 
below  the  electrode  tip  position,  between  the  soma  and  the  streak.  The  stippling  effect 
in  the  image  is  from  ganglion  cell  and  perhaps  displaced  amacrine  cell  somata  . 
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Appendix  B 

Whole- Cell  Patch 
Recordings  in  Isolated 
Intact  Turtle  Retina 


In  order  to  investigate  the  synaptic  inputs  to  DS  retinal  ganglion  cells  we 
have  developed  a  whole-cell  patch  technique  in  intact  isolated  turtle  retina. 
This  method  allows  the  following: 

•  Initial  extracellular  deterroination  of  receptive  field  properties,  with 
on-cell  patch  recordings. 

•  Intracellular  recordings  that  are  stable  for  typically  >  2  hours,  and 
that  are  quieter  than  conventional  microelectrode  recordings  due  to 
the  low  impedance  of  the  whole-cell  electrodes. 

•  Dialysis  of  the  cell  interior  by  the  electrode  solution  allowing  manipu¬ 
lation  of  the  intracellular  environment. 

In  this  appendix  we  shall  present  the  protocol  used  and  basic  properties 
of  cells  measured  with  the  whole-cell  patch. 

B.l  Recording  Method 

We  shall  now  describe  the  whole-cell  patch  protocol  materials  and  methods. 
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B.1.1  Electrodes  and  Electrode  Solutions 

Whole- cell  patch  electrodes  are  made  from  Coming  7052  capillary  glass, 
1.65  nim  O.  D.  1.20  mm  I.  D.  (A-M  Systems,  Everett,  Washington)  using 
a  three  pull  sequence  on  a  Flaming-Brown  micropipette  puller  (Model  P- 
80/PC,  Sutter  Instrument  Corporation,  San  Rafael,  California).  A  typical 
sequence  is  as  follows: 


Cycle 

Heat 

Pull 

Velocity 

Time 

1 

400 

0 

30 

80 

2 

400 

0 

30 

80 

3 

410 

60 

30 

40 

The  electrodes  have  a  steeply- tapering  tip  (diameter  approximately  2 
nm)  and  a  resistance  of  approximately  4  Md  when  filled  with  the  solutions 
to  be  described.  We  have  not  found  it  necessary  to  coat  the  electrodes  nor 
otherwise  treat  them  further  other  than  a  simple  fire-polishing  using  a  single 
(extra)  heat  cycle  of  the  micropipette  pullei.  This  step  consists  of  moving 
the  tip  of  one  electrode  at  a  time  into  the  neating  element  area,  while  the 
electrode  is  still  mounted  in  the  puller.  The  puller  is  then  activated  (Heat  = 
350,  Time  =  100)  while  holding  the  tip  in  this  position.  A  typical  electrode 
is  shown  in  Figure  B.l. 

Intracelluar  solutions  are  derived  from  the  protocol  of  the  first  whole¬ 
cell  patch  recordings  in  intact  vertebrate  retina  reported  by  Coleman  and 
Miller  [19].  Control  solutions  are  composed  of  90  mM  KCH3SO4,  5  mM 
NaCH3S04,  12  mM  MgCU,  1  mM  CaCU,  11  mM  EGTA,  5  mM  HEPES, 
2  mM  glucose,  and  1  mM  ATP,  titrated  with  KOH  to  a  pH  of  7.4.  For 
protocols  which  include  the  presiuned  blocking  of  GABAa  responses  of  the 
recorded  cell  alone  the  control  solution  is  used  minus  the  ATP  ([60])  and, 
occasionedly,  the  MgCh.  All  the  results  presented  in  this  thesis  are  from 
protocols  using  the  control  solution  without  ATP. 

B.l. 2  Electrode  Manipulation 

Electrodes  are  backfilled  with  a  hypodermic  needle  and  syringe,  which  is 
htted  with  a  2/im  in-line  filter.  Following  the  recommendation  in  [19],  the 
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Figure  B.l:  Whole-cell  patch  electrode  made  with  pull  sequence  in  text  (magnification 
570x).  Rtieetrodt  =  4Mn  when  filled  with  solution  described  in  text.  Tip  diameter 
«  2/im. 


hypodermic  needle  is  flushed  with  solution  prior  to  filling  an  electrode.  The 
filled  electrode  is  then  placed  in  a  patch  electrode  holder  (Model  PC-S3, 
E.W.  Wright,  Guilford,  Connecticut)  which  had  a  side  port  for  pressure 
communication  with  the  inside  of  the  electrode.  Tygon  tubing  (3/16”  0. 
D.  3’  length)  coimects  the  electrode  holder  to  a  lOcc  syringe  via  a  three- 
way  valve,  with  one  leg  of  the  valve  open  to  the  atmosphere.  This  syringe 
is  used  to  apply  either  pressure  or  vacuum  to  the  inside  of  the  electrode. 
The  holder  is  mounted  on  a  electrical  micropositioner  drive  (Nanostepper 
Model,  Scientific  Precision  Instruments,  Neckargemund,  FRG).  0^:^  the 
electrode  has  been  mounted  in  the  holder,  strong  pressure  (approximately 
7cc  of  movement  on  the  syringe)  is  applied.  This  pressure  has  three  pre¬ 
sumed  fimctions;  first,  to  force  intracellular  solution  to  the  tip  in  order  to 
finish  filling  the  electrode,  second,  to  cause  an  outflow  to  avoid  tip  clog¬ 
ging  as  the  tip  enters  the  bath  solution,  and  third,  to  prevent  the  tip  from 
getting  fouled  by  either  vitreous  or  extracellular  components  prior  to  being 
juxtaposed  against  a  cell  membrane. 

After  pressme  is  applied,  the  three-way  valve  is  positioned  to  seal  off  the 
electrode  from  the  other  two  outlets,  and  the  electrode  is  rapidly  advanced 
into  the  bath  solution.  The  amplifier  (Dagan  Model  3900,  Dagan  Corpora¬ 
tion,  Minneapolis,  Minnesota),  is  set  to  current  clamp  mode  and  adjusted 
to  compensate  for  the  the  electrode  series  resistance  sind  offset  potential. 
Blocked  electrodes  or  those  whose  resistance  was  not  between  3-6Mn  are 
discarded.  The  electrode  tip  is  then  lowered  onto  the  retina  under  video 
observation,  with  a  1  nA  40  Hz  square-wave  applied  current.  The  angle  be¬ 
tween  the  electrode  and  the  retina  surface  is  %45  degrees.  The  electrode  is 
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lowered  in  the  Nanostepper  continuous  mode  until  contact  with  the  vitreous 
is  detected  visually. 

B.1.3  Electrode  Solution  Bleeding 

After  contact,  electrode  advancement  is  continued  in  the  Nanostepper’s  “De¬ 
fined  Step”  mode  (3^m  step,  maximum  adjustment  of  speed  and  current 
potentiometers,  speed  switch  at  lOx,  1ms  delay)  until  either  the  electrode 
resistance  increases  suddenly  or,  more  t}rpically,  a  sudden  “bleeding”  of  the 
tip  solution  into  the  tissue  is  observed  on  the  monitor.  At  first  the  dramatic 
appearance  of  a  large  (roughly  100  -  200  fan  diameter)  stain  was  thought  to 
be  due  to  tip  breakage,  but  subsequent  inspection  of  the  electrode  showed 
that  this  was  not  the  case.  Apparently  the  bleeding  occurs  when  the  elec¬ 
trode  tip  leaves  the  extra-retinal  space  (including  the  dense  vitreous)  and 
enters  the  ganglion  cell  layer  (GCL).  It  is  possible  that  the  resistance  to 
flow  of  solution  out  of  the  electrode  tip  is  reduced  in  the  GCL,  and  that  flow 
is  also  constrmned  within  a  very  narrow  layer.  This  may  explain  why  an 
apparently  very  small  amount  of  electrode  solution  could  make  such  a  large 
visual  impression.  In  any  event,  electrode  bleeding  has  correlated  strongly 
xmth  achieving  giga~seal  and  whole^cell  access. 

B.1.4  Making  the  Seal 

Immediately  when  the  bleeding  occurs  the  three-way  valve  is  quickly  moved 
to  open  communication  between  the  electrode  and  the  syringe.  The  goal 
here  is  to  reduce  the  positive  pressure  significantly  within  the  electrode  (as 
evidenced  by  the  subsequent  outward  movement  of  the  syringe)  without 
removing  the  pressure  entirely.  Reduction  of  pressure  is  often  accompanied 
by  an  increase  in  the  tip  resistance,  typically  several  Mfl.  Occasionally  the 
giga-seal  forms  spontaneously  at  this  point.  If  not,  the  tip  is  slowly  advanced 
(same  Nanostepper  settings  as  above,  but  speed  switch  set  to  lx)  until  the 
tip  resistance  increases  to  approximately  10  -  20  MO.  The  three-way  valve 
is  then  switched  so  that  the  inside  of  the  electrode  is  exposed  to  atmospheric 
pressure.  If  a  giga-seal  does  not  form  at  this  point  on  its  own,  slight  suction 
is  applied  either  by  mouth  or  by  switching  the  three-way  valve  to  re-connect 
the  syringe  with  the  electrode  and  moving  the  syringe  very  slightly  (<  <lcc), 
until  the  gigsi-seal  starts  to  form.  When  the  giga-seal  forms,  as  evidenced  by 
the  rapid  increase  in  the  tip  resistance  (tip  voltage  from  the  amplifier)  the 
current  clamp  waveform  is  qtiickly  shut  off  (0  current  clamp  mode).  This 
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avoids  large  voltage  excursions  which  can  hinder  seal  formation  or  cause 
premature  membrane  rupture. 

On-cell  patch  giga-seals  are  then  sdlowed  to  stabilize  for  typically  five  to 
ten  minutes  before  continuing.  Typically  seals  start  out  at  IGH,  and  in¬ 
crease  to  4  -  20Gn.  Seal  resistance  is  checked  periodically  by  100-200  msec 
10  pA  pulses  (see  Figures  B.2  and  B.3).  The  stain  from  the  earlier  bleeding 
typically  clears  completely  within  five  minutes.  Also,  the  apparent  ^Erett” 
in  the  on-cell  patch  mode  can  vary  from  -20  to  -|-5mV,  although  between  -5 
and  OmV  is  most  common.  It  is  not  known  what  causes  these  variations,  but 
possible  candidates  include  inconsistent  electrode  wire  (Ag/AgCl)  prepara¬ 
tion. 


B.1.5  On-Cell  Patch  Recordings  Prior  to  Whole-Cell  Access 

In  most  protocols  which  resulted  in  successful  whole-cell  recordings  of  spik¬ 
ing  cells,  high  signal  to  noise  ratio  extracellular  recordings  could  be  obtained 
when  the  on-cell  patch  ^ga-seal  stabilized  (Figure  B.2).  In  fact,  good  extra¬ 
cellular  recordings  are  almost  diagnostic  as  far  as  the  likelihood  for  successful 
whole-cell  access. 

When  on- cell  recordings  were  possible,  receptive  field  spatial  and  DS 
properties  were  evaluated  according  to  the  protocols  described  in  Section 
3.2.  As  mentioned  in  this  section,  the  success  of  these  protocols  was  variable 
due  to  the  var;;  ing  quality  and  stability  of  the  on-cell  patch.  Typical  on-cell 
recording  times  ranged  from  1  to  20  minutes. 

B.1.0  Getting  Whole-Cell  Access 

After  on-cell  patch  recordings  (if  any),  a  1  Hz  flashing  large  spot  stimulus  was 
turned  on  and  very  weak  suction  was  applied  by  mouth  or  with  the  syringe. 
Whole  cell  access  was  marked  by  either  a  fall  in  the  tip  potential  and/or 
photoresponses  (see  Figures  B.2  and  B.3).  If  the  cell  was  a  ganglion  cell 
spikes  were  often,  but  not  always,  seen  as  the  membrane  ruptured.  Stable 
resting  potentials  and  photoresponses  were  then  recorded  for  times  ranging 
from  30  minutes  to  over  4  hours. 

B.1.7  Importance  of  Visualization  of  Tissue  During  Record¬ 
ing 

In  our  experience  visualization  of  the  retina  and  electrode  is  vital  to  achiev¬ 
ing  whole-cell  recordings.  The  recordings  presented  here  were  done  at  low 
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Figure  B.2:  Transition  between  on-cell  patch  and  whole-cell  patch.  On-cell  patch  was 
monitored  by  hyperpolarizing  current  steps  (10  pA,  70  msec),  occurring  here  (top  panel) 
at  1.2,  2.4,  and  3.5  seconds.  In  the  top  panel,  whole-cell  access  was  indicated  by  drop 
in  Erttt,  light-evoked  spikes,  and  a  drop  in  Ri„,  Momentary  drops  in  Freiti  example 
as  seen  at  the  beginning  of  the  trace,  were  typical  when  suction  was  applied  to  the 
electrode.  After  whole-cell  access  is  achieved,  bursts  of  spikes  are  seen  in  response  to  a 
flashing  Ihz  spot  (ON/OFF  transitions  at  0,  500,  1000  ms,  etc.).  In  the  bottom  panel 
the  expansion  of  the  trace  in  top  panel  (starting  at  the  1  second  point  of  the  this  trace) 
shows  a  light-evoked  on-cell  patch  spike,  followed  by  the  current  step  response.  Rteai 
for  this  recording  is  about  5Gft,  and  Rjn  is  about  iGfl. 
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Figure  B.3;  The  on-cell  patch  to  whole-cell  patch  transition  in  this  cell  clearly  shows 
the  typically  4  to  10-fold  drop  in  Rjn  seen  with  the  transition  (10  pA,  70  msec  current 
steps).  A  flashing  light  stimulus  at  1  hz  throughout  this  recording  fails  to  generate 
extracellular  spikes,  which  is  confirmed  by  the  subthreshold  light-evoked  EPSPs  seen 
after  whole-cell  access.  The  burst  of  spikes  at  3400  ms  is  in  response  to  a  depolarizing 
current  step.  R,e<d  for  this  recording  is  about  8Gn  and  Rin  is  about  IGO. 
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magnification  (w  80x ),  where  individual  somata  were  just  at  the  limit  of  res¬ 
olution  (see  Figure  A.6).  Thus,  these  recordings  were  done  partially  blind, 
in  that  there  was  no  attempt  to  “stab”  a  specific  cell.  However,  completely 
blind  manipulation  of  the  electrode  (e.g.  tip  position  monitored  solely  by 
tip  resistance),  has  been  much  less  successful,  especially  since  the  bleeding 
described  earlier  cannot  be  followed. 

Also,  the  optics  of  our  setup  allow  high-magnification  ( >  500  x )  of  suf¬ 
ficient  quality  to  resolve  specific  somata  with  an  adjacent  electrode  tip  (see 
FigureB.4).  It  is  likely  that  with  minor  improvements  in  the  setup,  it  will 
be  possible  to  record  from  specific  cells,  e.g.  from  previously-labelled  cells. 


Figure  B.4:  Scanned  infra<red  video  hardcopy  image  of  intact  isolated  turtle  retina  at 
high  resolution  (magnification  »  800x)  made  during  an  experiment.  Cell  bodies  are 
clearly  visible  with  no  staining.  A  whole-cell  electrode  tip  is  seen  in  the  middle  of  the 
micrograph  approaching  a  soma.  The  shank  of  the  electrode  extends  from  the  upper 
left  corner  of  the  micrograph. 
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B.2  Whole*Cell  Patch  Recordings  Show  Both  Light 
and  Current-Evoked  Spikes,  Long  Time  Con¬ 
stants,  and  High  Input  Resistances 

The  whole- cell  patch  recordings  show  stable  and  vigorous  action  potentials 
in  response  to  both  flashing  light  stimuli  and  depolarizing  current  pulses, 
as  shown  in  Figure  B.5.  This  figure  also  demonstrates  the  typical  long  time 
constimt  (r  for  the  dashed  current  step  response  above  is  on  the  order  of  50 
milliseconds)  and  high  input  resistance  (700  MH  for  this  cell)  of  cells  using 
the  whole-cell  technique. 


Figure  B.5:  Whole-cell  patch  recordings  of  an  ON/OFF  ganglion  cell,  showing  responses 
to  both  flashing  spot  stimuli  and  depolarizing  current  pulses.  On/off/on/off  transitions 
of  light  spot  are  at  0,  500,  1000  and  1500  ms,  respectively.  On/ofF/on/ofF  transitions 
of  current  pulse  are  at  -100,  750,  1150  and  1900  ms,  respectively.  Non-uniform  spike 
heights  are  partly  due  to  the  sampling  rate  of  1/ms. 


We  believe  that  the  recordings  presented  here  are  likely  to  be  from  gan¬ 
glion  cells  for  two  reasons:  1)  the  vigorous  spike  responses  are  more  typical 
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of  the  classic  gcinglion  cell  response,  rather  th2Ui  from  amacrine  cells,  and 
2)  the  electrode  seal  is  established  very  soon  after  initial  contact  with  the 
retinal  surface  (ganglion  cell  layer).  However,  it  is  still  possible  that  some 
recordings  are  from  spiking  displaced  2imacrine  cells. 


B.3  Reduction  of  IPSPS  Over  Time  With  ATP- 
Less  Whole-Cell  Patch  Electrode  Solution 

In  most  of  the  retinal  cells  that  we  have  recorded  from  using  the  whole-cell 
technique,  IPSPs  are  observed  for  large  spots  (typically  >  200/im)  for  the 
first  10  to  30  minutes  of  the  recordings.  When  ATP  is  not  included  in  the 
electrode  solution  the  IPSPs  eventually  disappear,  often  losing  more  than 
50%  of  amplitude  within  a  30  second  period.  EPSPs,  on  the  other  hand, 
are  much  more  stable,  with  many  cells  showing  less  than  10%  attenuation 
of  EPSPs  over  1  to  2  hours.  Figures  B.6  through  B.12  demonstrate  the 
reduction  of  IPSPs  to  large  spots  over  time. 

B.3.1  Antagonistic  Center/Surround  is  Maintained  Despite 
Local  Reduction  of  IPSPS 

That  the  attenuation  of  inhibition  using  our  technique  is  a  loco/ phenomenon 
is  supported  by  the  maintenance  of  an  inhibitory  siirround  in  most  cell^. 
This  is  demonstrated  in  Figures  B.IO,  B.ll,  and  B.12,  where  the  largest 
response  continues  to  be  from  a  sub-maximal  spot,  despite  lack  of  IPSPs. 
These  results  also  point  out  that  antagonistic  center/surround  properties 
are  derived  from  interactions  both  pre  and  post-synaptic  to  the  ganglion 
cell,  and  is,  in  some  cases,  purely  post-synaptic. 


B.4  Measurement  of  Passive  Parameters  with  Whole 
Cell  Patch  Recordings 

Data  from  59  cells  recorded  under  whole-cell  conditions  were  pooled  to  es¬ 
timate  Rjn,  To,  and  Ere«(-  The  protocol  used  was  a  10  trial  average  of  the 
response  under  current  clamp  to  a  10  pA  square  wave  (no  DC  offset).  Lin¬ 
ear  regression  on  the  logarithm  of  this  response  was  used  to  derive  tq.  An 
example  response  is  shown  in  Figure  B.13. 
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Figure  B.6:  Light  responses  of  ON  ganglion  cell  (A29)  to  1200  fim  diameter  spot  at 
times  noted,  relative  to  start  of  whole-cell  access.  These  traces  are  averages  of  4  trials. 
On/off/on/off  transitions  of  light  spot  are  at  0,  500,  1000  and  1500  ms,  respectively. 
Recording  at  17  minutes  is  subthreshold  for  spikes,  and  spikes  were  lost  in  this  cell  after 
30  minutes.  In  this  particular  cell,  the  resting  potential  depolarized  over  time,  but  it 
was  more  typical  to  sec  a  drop  in  Ertu  during  the  recording.  IPSPs  relative  to  Ere$t 
are  seen  at  light  OFF  in  the  earlier  recording,  but  this  negative  phase  is  lost  in  the  later 
recording.  Note  that  the  rise  in  £re«t  thi*  cell  would  be  expected  to  amplify  any 
IPSP. 
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Figlire  B.7:  Responses  of  ganglion  cell  A84  to  400  and  800  nm  diameter  spot  at  start  of 
whole-cell  access  (average  of  4  trials).  On/off/on/off  transitions  of  light  spot  are  at  0, 
500,  1000  and  1500  ms,  respectively.  Stimuli  are  subthreshold  for  spikes.  Here  the  cell 
shows  a  center/surround  OFF  response,  with  a  complicated  (negative)  ON  response. 
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Figure  B.8:  Response  of  cell  AJ» taken  more  than  one  hour  after  start  of  whole-cell 
access  (average  of  4  trials).  Er„t  =  -55  mV.  This  cell  now  show  an  ON/OFF  response, 
and  each  edge  response  is  simpler  than  the  initial  recording.  These  responses  suggests 
that  inhibitory  input  to  the  cell  has  been  blocked,  uncovering  an  underlying  excitatory 
ON/OFF  input.  Actual  blockage  of  inhibitory  synaptic  channels  is  further  supported  by 
the  slight  rise  in  Ert»t  5mV),  although,  as  mentioned  earlier,  it  was  more  typical  to  see 
a  drop  in  Erttt  over  time.  Note  also  the  loss  of  antagonistic  center/surround  properties 
in  this  particular  cell.  As  noted  in  the  text  and  in  Figures  B.IO  to  B.12,  in  many  cells 
a  loss  of  IPSPs  was  seen  over  time  with  nonetheless  maintenance  of  center/surround 
properties.  As  before,  on/off/on/off  transitions  of  light  spot  are  at  0,  500,  1000  and 
1500  ms,  respectively. 
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Figure  B.9:  Response  of  cell  A84  taken  immediately  after  that  in  Figure  B.8,  with 
addition  of  depolarizing  current  (to  -40  mV)  in  order  to  unmask  any  “silent"  IPSPs 
(average  of  4  trials).  No  negative  phase  is  observed,  confirming  the  blockage  of  inhibitory 
input  over  time. 
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Figure  B.IO:  ON  and  OFF  responses  of  ganglion  cell  AlOO,  3  minutes  after  star^  of 
whole-cell  access,  to  1  hz  50  to  1600  ftm  diameter  spots  (average  of  2  trials).  Ere$t  = 
-59  mV.  On  and  off  transitions  of  light  arc  at  0  and  500  ms,  respectively.  This  cell  shows 
a  strong  OFF  component,  with  a  weaker  ON  component,  each  of  which  is  strongest  for 
the  100-200  /im  spots.  An  IPSP  to  both  ON  and  OFF  transitions  is  seen  for  the  larger 
spots. 
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Figure  B.ll:  Responses  of  ceil  AlOO  under  same  conditions  as  Figure  B.IO,  31  minutes 
after  start  of  whole<celi  access.  The  antagonistic  center/surround  behavior  is  main¬ 
tained,  but  only  slight  IPSPs  to  the  larger  spots  are  seen.  However,  the  resting  potential 
has  dropped  by  about  10  mV  {Er„t  =  mV),  which  may  hide  the  relatively  small 
IPSPs  seen  in  Figure  B.IO. 
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Figure  B.12:  Responses  of  cell  AlOO  under  same  conditions  as  the  previous  two  fig¬ 
ures,  but  with  a  depolarizing  holding  current  (to  -58  mV)  in  order  to  “unmask”  any 
hidden  IPSPs  of  Figure  B.ll.  However,  the  IPSP  phase  seen  in  Figure  B.IO  is  not 
observed  here,  suggesting  that  IPSPs  have  been  blocked.  The  fact  that  an  antagonis¬ 
tic  center/surround  characteristic  remains  suggests  that  the  inhibition  of  the  inhibition 
is  local  to  the  recorded  cell.  Note  the  spike-induced  fast-afterhyperpolarization  in  the 
OFF  response  to  the  200^m  spot,  suggesting  that  voltage-dependent  re/hyperpolarizing 
mechanisms  are  preserved. 
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Figure  B.13:  Example  current  clamp  record  of  response  to  10  pA  square  wave  (average 
of  10  trials).  Erttt  =  *S5  mV.  The  results  in  this  section  are  derived  from  similar 
responses  of  59  cells. 


B.4.1  Estimates  of  iZ/„  and  tq 

We  h<ive  lound  that  the  whole-cell  patch  technique  with  turtle  retinal  gan¬ 
glion  cells  gives  relatively  high  values  of  Rin  (mean  =  O.STGH)  and  tq  (mean  = 

56  milliseconds).  The  distributions  of  these  parameters  are  shown  in  Figures  B.14 
and  B.I5. 

These  values  are  consistent  with  reports  in  other  neurons  using  the 
whole-cell  patch  technique.  For  example,  Coleman  rmd  Miller  ([19])  give 
mean  values  of  iZ/n  =  l.SGH  and  tq  =  €8  ms.  Likewise,  in  the  salamander 
retinal  slice  preparation,  Lukasiewicz  and  Werblin  ([46])  report  mean  values 
of  Rin  =  2.9Gn  and  tq  =  41.3  ms.  Similar  values  have  been  reported  from 
non-retinal  netirons  as  well. 

If  we  assume  that  the  specific  capacitance  Cm  =  1.0/iF/cm*,  then  the 
measured  tq  of  56  milliseconds  implies  a  specific  resistivity  of  Rm  =  56Kncm^. 
This  value  is  consistent  with  estimates  from  at  least  two  modelling  studies 
(cerebellar  Purkinje  cells  [58],  Rm  >  45Kflcm’;  hippocampal  pyramidal 
cells  [9],  Rm  ^  dOKHcm^)  which  derived  Rm  by  various  and  independent 
methods. 
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Figure  B.14:  Distribution  of  Rjn  values  for  a  series  of  59  cells.  The  mean  value  for 
Rin  =  0.87Gn  (SEM  =  0.09). 
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Figure  B.15:  Distribution  of  ro  values  for  the  same  59  cells  shown  in  Figure  B.14.  The 
mean  value  of  t))  =  56  milliseconds  (SEM  =  3). 
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B.4.2  Effect  of  Non-Linear  Membrane  on  i2/n  and  tq 

Given  the  fairly  large  voltage  excursion  (typiceJly  10  mV)  and  long  time 
constant  of  the  step  responses,  and  that  the  responses  had  both  a  hjrper- 
polarizing  and  depolarizing  component  (relative  to  the  resting  potential),  it 
is  likely  that  non-linear  conductances  participated  in  shaping  the  responses. 
Since  many  rectifying  channels  tend  to  activate  in  this  voltage  range  (partic¬ 
ularly  for  depolarizations),  it  is  possible  that  this  protocol  underestimated 
both  Rin  and  tq. 

B.4.3  Effect  of  Electrode  Seal  on  H/n  and  tq 

In  addition,  given  the  the  higi.  values  for  both  and  tq  it  is  useful  to 
consider  the  contribution  of  the  giga-seal  (between  the  electrode  and  the 
extracellular  space)  on  these  parameters.  A  rough  estimate  would  be  to 
assume  a  measured  value  of  Rjn  =  IGH  and  consider  a  lower  bound  of 
^itai  =  2Gn.  For  this  value  of  Rttai^  balf  of  the  input  conductance  (l/Rin) 
would  be  due  to  electrode  leak,  and  thus  the  true  jR/„  =  2GO.  Likewise,  To 
under  these  conditions  would  be  underestimated  by  a  factor  of  2.  In  practice, 
Rteal  >  SGH,  but  this  value  is  of  course  measured  prior  to  whole-cell  access, 
and  it  is  possible  that  the  seal  is  partly  compromised  by  the  whole-cell  access. 
In  summary,  these  factors  may  also  contribute  to  an  underestimation  of  both 
Rin  and  tq. 


B.4.4  Estimation  of  Erett 

After  IPSP  attenuation  (see  Section  B.4.5),  the  resting  potential  tended  to 
be  stable  for  the  remainder  of  the  recording.  For  the  set  of  59  cells  just 
presented,  the  mean  value  for  Erett  was  -56  mV.  The  distribution  of  Erett 
values  in  these  cells  is  shown  in  Figure  B.16. 

B.4.5  Correlations  Between  Passive  Parameters 

As  shown  in  Figures  B.17,  B.18  and  B.19,  there  was  only  a  slight  positive 
pair-wise  correlation  between  Rjnt  fo  and  Ereto  One  possible  mechanism 
contributing  to  these  correlations  involves  the  state  of  the  inhibitory  synap¬ 
tic  channels.  Typically,  the  attenuation  of  IPSPs  during  whole-cell  access 
was  accompanied  by  a  drop  in  Erett  by  «  10  to  15  mV.  In  a  few  cells, 
Rjn  also  dropped  as  the  IPSPs  disappeared  by  approximately  50%.  These 
preliminary  observations  suggest  that  blockage  of  IPSPs  may  be  caused  by 
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tonic  activation  of  the  associated  synaptic  channels,  when  there  is  no  ATP 
in  the  electrode  solution. 


B.5  Staining  of  Recorded  Cells 

We  have  also  done  trials  of  whole-cell  recording  with  Lucifex  Yellow  in  the 
electrode  solution  (not  shown).  Preliminary  results  indicate  that  this  stain¬ 
ing  method  works  well  in  our  turtle  preparation. 
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Figure  B.16:  Distribution  of  values  for  the  same  59  cells  shown  in  Figure  B.14. 
The  mean  value  of  Ere$t  =  —56  mV  (SEM  =  1). 


Figure  B.17:  Scatter  plot  of  Rin  and  tq  values  for  the  59  cells  shown  in  Figure  B.14. 
Solid  line  is  least-squares  linear  fit  (correlation  coefficient  r  =  0.13). 
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Figure  B.18:  Scatter  plot  of  iZ/n  and  Ertti  values  for  the  59  cells  shown  in  Figure  B.14. 
Solid  line  is  least-squares  linear  fit  (correlation  coefficient  r  =  0.33). 


Figure  B.19:  Scatter  plot  of  and  Erft  values  for  the  59  cells  shown  in  Figure  B.14. 
Solid  line  is  least-squares  linear  fit  (correlation  coefficient  r  =  0.14). 
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Appendix  C 

The  SURF-HIPPO  Neuron 
Simulator 


C.l  Introduction 

The  SURF-HIPPO  Neuron  Simulator  is  a  circuit  simulation  package  for 
investigating  morphometricly  and  biophysically  detailed  models  of  single 
neurons  and  small  networks  of  neurons  (Figure  C.l).  SURF-HIPPO  is  based 
in  part  on  the  SURF  circuit  simulator,  written  by  Don  Webber,  VLSI  CAD 
Group  at  the  Un  versity  of  CiUfomia  at  Berkeley.  SURF  was  written  in 
both  Common  LISP  and  "'LISP,  and  SURF-HIPPO  retains  this  option  (the 
'"LISP  version  for  the  Connection  Machine  (CM2)  is  under  development). 

SURF-HIPPO  allows  ready  construction  of  multiple  cells  from  various 
file  formats  (including  the  format  used  by  the  simulator  package  NEURON, 
written  by  Hines,  see  below),  which  can  describe  complicated  dendritic  trees 
in  3-space  with  distributed  non-linearities  and  synaptic  contacts  between 
cells.  The  retention  of  XYZ  coordinates  for  each  circuit  node  may  be  ex¬ 
ploited,  for  example,  by  simulations  which  use  spatially- coded  input,  or  with 
multiple  cell  simulations  that  have  contacts  between  cells  defined  by  spatial 
proximity  (Figure  C.2).  Stationary  and  moving  two  dimension^d  light  input 
(for  retinal  simulations)  is  also  provided.  Multi- compartment  intracellular 
calcium  accumulation  is  provided  for  modelling  calcium- dependent  processes 
(e.g.  C'a'*"'‘-dependent  channels). 

This  chapter  provides  an  overview  of  SURF-HIPPO.  Complete  docu¬ 
mentation  is  included  with  the  program  code. 
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C.2  Numerical  Methods 


The  serial  version  of  SURF-HEPPO  uses  a  numerical  technique  adapted  from 
that  of  Hines  ([37])  with  some  changes,  including  adaptive  time  steps.  *LISP 
simulations  use  the  relaxation-based  methods  exclusively.  At  present,  all 
simulations  are  done  with  the  direct  (serial)  methods. 

Briefly,  the  key  contribution  of  the  Hines  algorithm  is  an  elegant  ordering 
of  tree  circuit  topologies  -  clearly  relevant  for  neurons  -  that  results  in  a  tri¬ 
diagonal  circuit  matrix.  The  resulting  matrix  equation  may  then  be  solved 
in  0{n)  time,  as  opposed  to  the  O(n^)  time  for  a  general  matrix.  Hines 
also  proposed  that  the  solving  of  voltage-dependent  Hodgkin-Huxley  channel 
particle  states  may  be  done  out  of  the  loop  (i.e.  at  times  interleaved  between 
the  node  voltage  time  steps),  since  there  is  an  inherent  time  delay  between  a 
given  node  voltage  and  the  state  of  any  channel  particles  controlled  by  that 
node.  This  explicit  solution  of  the  particle  states  simplifles  and  speeds  up 
the  overall  circuit  simulation. 


C.3  Channel  Models 

All  voltage-dependent  channels  are  described  with  a  modified  Hodgkin- 
Hiixley  model  ([10]).  An  advantage  of  this  formulation  is  that  channel  par¬ 
ticle  parameters  are  consister.t  and  may  be  related  to  specific  biophysical 
mechanisms. 


C.4  System  Versions 

One  version  of  SURF- HIPPO,  complete  with  extensive  user  interface  (menus, 
3-D  graphics  of  dendritic  trees,  automatic  plotting)  runs  in  the  Symbolics 
LISP  Machine  environment.  Another  (with  a  user  interface  currently  tmder 
development)  nms  in  the  Sun  environment  (under  Austin  Kyoto  Common 
Lisp  (AKCL)).  SURF-HIPPO  is  not  ’’public  domain”  in  the  sense  that  it  is 
not  supported,  but  inquiries  are  welcome. 


C.5  SURF-HIPPO  Basics 

Functions  which  describe  circuits  may  be  written  in  a  variety  of  different 
ways,  allowing  for  a  flexible  and  relatively  general  input  base  for  circuit 
descriptions.  The  simplest  method  is  by  calling  the  appropriate  sequence 
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of  “create-segment”  functions.  A  more  sophisticated  method  of  construct¬ 
ing  (complicated)  cell  geometries  is  with  the  ”create-tree”  function,  which 
takes  a  list  of  tree  node  parameters  and  generates  the  "create-segment”  calls 
automatically. 

Another  method  of  entering  cell  descriptions  is  by  tracing  the  histol¬ 
ogy  directly  on  the  screen,  using  SURF-HEPPO’s  Dynamic  Neuron  Builder 
(Figure  C.3). 

Finally,  as  mentioned  above,  SURF-flIPPO  can  read  cell  description  files 
in  the  format  used  for  Hines’s  NEURON  program. 

C.6  Representation  of  Circuit  Structure 

Symbolically,  circuits  are  composed  of  a  set  of  object  instances,  of  which 
there  are  several  different  types:  nodes,  segments,  somata,  channels,  synapses, 
current  sources,  particles,  concentration  integrators,  etc.  There  are  many 
different  interrelationships  between  objects  in  a  hierarchical  fashion,  e.g.  an 
instance  of  a  channel  object  references  an  instance  of  a  node  object  (corre¬ 
sponding  to  the  electrical  circuit  node  that  it  is  attached  to),  an  instance  of 
a  segment  or  soma  object  that  it  is  a  part  of,  instances  of  particle  and/or 
concentration  particle  objects  that  determine  the  channel's  conductance, 
and,  indirectly,  instances  of  particle  nodes  which  correspond  to  the  parti¬ 
cles’  states. 

The  structure  organization  for  a  typical  circuit  element  is  as  follows: 

;S«gB«nt  nan*  it  saB*  a  distal  nod#  naaa.  Segments  are  added  from  soma  out. 
(def struct  (segment  :conc-name) 

(name  ""  ’.type  string) 

g-axial 

;  Proximal  node  Prox  o — /\/\/\ - ♦ — o  Distal 

; Distal  node  i 

:  memb-elements 

:  I 

:  Gnd 

(node-1  nil  ’.type  node) 

(node-2  nil  :type  node) 

(cell  nil  ’.type  cell) 

;For  retinal  cells:  1  to  6  for  dendrites  in  the  inner 
;plexiform  layer.  IIL  if  not  in  IPL. 

(ipl-stratum  nil  s'*:/?*  integer) 
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: microns 

(length  zero  .‘type  long-float) 

(diameter  zero  :type  long-float) 

; Elevation  angle  relative  to  proximal  branch  orientation. 

(theta  zero  long-float) 

: Azimuth  angle  relative  to  proximal  branch  orientation. 

(phi  zero  ^'type  long-float) 

(core-seg  nil  :type  core-segment) 

(model  nil  :type  model- instance) 

:  microsiemens 

(g-axial  zero  :type  long-float) 

(g-leak  zero  :type  long-float) 

;  millivolts 

(v-leak  -70.0  ‘.type  long-float) 

; nanofarads 

(capacitance  zero  :type  long-float) 

; Concentration  integrator 
(conc-int  nil) 

;Zf  this  segment  is  connected  to  a  Hines  node,  then  this  is  the 
{index  for  that  node. 

(branch-node-index  nil)) 

Simulation  of  circuits  is  in  terms  of  their  nodes,  which  include  the  actual 
electrical  nodes  (corresponding  to  somata  and  segments)  as  well  as  elements 
whose  states  are  both  dependent  on  other  nodes  and  in  turn  control  either 
br^mch  elements  of  the  circuit  or  other  (non- voltage)  nodes.  This  latter  class 
of  elements  include  the  following: 


Element  Type 

"Voltage”  = 

Dependent  On 

Controls 

V-dep  particles 

Prob.  of  open  state 

Membrane  voltage 

Sehannel 

Ca+'^'-dep  particles 

Prob.  of  open  state 

Ca'^'^\  Integrator 

9channel 

ICa 

Ca'^'^'-dep  particles 

The  “voltages”  of  the  circuit  nodes  correspond  to  the  state  variables  of 
the  circuit.  Note  that  branch  elements  (resistors,  channels,  capacitances, 
sjmapses)  do  not  represent  the  state  variables  per  se;  their  values  are  com¬ 
puted  in  terms  of  the  appropriate  node  ”voltages”. 
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Figure  C.l:  An  example  SURF*HIPPO  simulation.  Various  windovi^s  on  the  screen 
show  node  voltages,  membrane  currents,  channel  particle  states,  and  cell  histology  with 
stimulus. 


Figure  C.2:  Cell  descriptions  entered  into  SURF-HIPPO  include  3-dimensional  node 
coordinates,  useful  for  spatially-coded  input  or  for  establishing  spatially-dependent  re¬ 
lationships  between  multiple  cells.  In  this  figure  the  simulated  star  amacrine  cell  from 
Figure  C.l  is  shown  at  a  viewing  angle  of  2°  relative  to  the  retinal  plane. 
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Figure  C.3:  The  Dynamic  Neuron  Builder  allows  the  entering  of  a  cell  geometry  directly 
on  the  screen,  using  the  mouse.  For  example,  in  this  figure  the  histology  of  a  putative 
OS  turtle  retina  ganglion  cell  ([38])  from  Kolb  (type  G20  [42])  was  xeroxed  onto  a 
transparency,  which  was  then  taped  to  the  computer  screen  and  traced  with  the  mouse. 
The  Dynamic  Neuron  Builder  also  allows  the  inclusion  of  various  synaptic  arrangements, 
specified  by  menu  commands.  In  this  figure,  excitatory  synapses  (filled  circles)  were 
automatically  added  in  an  annulus  centered  on  the  soma  (160  ftm  I.D.,  400/im  O.O.). 
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